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Preface 


For a long time, quality has been one of industry’s main preoccupations. 
It remains so today. 

There is some foundation for the statement that there is a ‘quality 
crisis’ in Europe, the methods traditionally used in our industries being 
unable to meet today’s demands. Consequently it is essential to look for 
new directions in which to progress, taking account of the methods for 
achieving quality that have been developed in recent years. These 
methods impact on all parts of the industrial enterprise — marketing, 
manufacturing, research and development, after-sales services. All staff, 
administrative or technical, are involved. 

The present book describes the tools that can help anyone who is 
concerned with the concept of ‘total quality’; it will also be a valuable 
educational aid for students reading for degrees or other qualifications 
in engineering. 


P. Lyonnet 
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General questions and 
concepts 


1 INTRODUCTION: HOW CAN WE ACHIEVE TOTAL 
OUALITY? 


This book is concerned with the various techniques and methods of 
analysis that can be used to ensure total quality in a project. In this first 
chapter we show the costs that result from not achieving quality, so as 
to make clear how important a quality-assurance service is to any 
enterprise. We stress also the involvement of marketing, particularly in 
laying down specifications for reliability. . 

Current techniques concerning reliability are developed in Chapter 2; 
here, much space is given to quantitative analysis, for this enables the 
whole range of problems that can be raised by questions of reliability to 
be dealt with in a co-ordinated manner. The advantage of using 
reliability as a guide to the choice of a technology is stressed. 

Chapter 3 is devoted to methods for controlling manufacturing proces- 
ses, in particular the implementation of control charts. 

Chapter 4, effectively a continuation of Chapter 3, deals with the 
management of input of raw materials and other supplies and output of 
finished products. All the methods proposed are described in sufficient 
detail to make immediate application possible. 

Chapter 5 develops methods for causal analysis — inferring causes 
from observed effects. The role of these methods in various parts of the 
enterprise is discussed, and their use in ‘quality circles’. 

Chapter 6 gives the mathematical apparatus needed to support the 
techniques described in the previous chapters, and the book concludes 
with a set of exercises, followed by their solutions; we trust that this will 
add to its educational value. 


L2" WHY QUATITY? 


For a long time, quality specialists were accused of being more con- 
cerned with formalities than with productive work; this is no longer the 
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case, but nevertheless it is important to emphasize areas where quality 
can prove a source of profit. 


1.2.1 Measures of ‘non-quality’ 


(a) Customer action 
If a product fails to meet its specification during the guarantee period, 
costs are incurred that are easily measurable: 


@ involvement of the after-sales service; 

@ costs of repair or modification; 

@ transport costs; 

® time wasted; 

@ possibility of penalty payments to the customer. 


The cost of not achieving quality can be determined from these. 


(b) Losses within the enterprise 

Losses within the enterprise are mainly the costs of scrapping or 
reworking the product; they are easily determined from the costs of 
materials and labour. 


Example 1 

A sheet metal workshop produces laminations for rotors and stators of 
electric motors, with a value of 1 F per kilogram. A monthly loss of 2 
tonnes, because of poor quality, would cost £3000, which would justify 
employing several staff purely for quality control. 


(c) Quality level imposed by the customer 

There are some markets in which a supplier is not allowed to compete 
unless he has a quality-control system that ensures that a stated level is 
reached. Such markets are aerospace, defence, nuclear reactors, na- 
tional electric power stations etc.; they can insist on 


@ certification by the appropriate standards body; 

© quality-assurance manuals, drawn up by a specialist organization; 

e manufacturing control documentation, to show that the required 
standards have been adhered to; 

@ development of a checking procedure. 


(d) Influence of quality on sales; index of quality 

Obviously the quality of a product will affect the demand for it. A 
customer who has been let down by the performance of a product is not 
a good advertisement for the manufacturer, and sales are likely to suffer 
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in consequence. Whilst it may not be easy to quantify this effect it is 
nevertheless of real importance, and some indicators have been devised 
to help justify the cost of a quality-assurance service. 

A conventional scale for denoting quality is shown in Fig. 1.1. The 
aim is to achieve an overall level of around 1 for a product. The best 
procedure is to estimate values of this index for the important features 
such as reliability, maintainability, aesthetics etc., with weighting factors 
indicating their relative importance, and to calculate an overall index 
from these; this will reveal the main weaknesses. Table 1.1 is an 
example. 


High quality 2 
Optimum quality 1 
Critical level 0.5 


Very poor quality 0 


Figure 1.1 Conventional scale for quality. 


Table 1.1 Weightings for separate features 


Estimated Weight Overall 
Feature index (%) index 
Reliability Led 100 
Maintainability 0.8 “ 80 
Performance 2; 100 
Support costs 0.9 90 
Aesthetics ile) 60 1.08 


(e) Quality-price-demand relations 
As we have noted, good quality will increase the demand for a product, 
poor quality will depress this; but demand will be depressed also by a 
high price. If a product sells badly because of poor quality the price will 
have to be lowered in order to increase sales. The situation can be 
represented by a set of curves as in Fig. 1.2. 

A study of these relations between quality, price and demand enables 
the importance of the quality of a product to be assessed. It is important 
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Volume of demand 


Figure 1.2 Quality versus price and demand. 


to recognize that for some products quality is the overriding considera- 
tion whilst for others the overriding consideration is price; thus one 
speaks of demand being inelastic with respect to quality or inelastic with 
respect to price. 


1.2.2 Evolution of the quality concept 


In the industrial world the concept of quality is a recent development; 
up to the end of the Second World War, in fact, it was scarcely taken 
into consideration. The various stages in its evolution often mirror the 
varying level of its adoption in industrial enterprises. We can distinguish 
the following chronological stages. 


(a) Production 

There was no quality service; everything was dominated by the manufac- 
turing process. To ensure that the products delivered conformed to their 
specifications a service independent of the production organization 
weeded out those that did not. 


(b) Statistical control 
Statistical control became common during the 1950s, particularly with 
the appearance of the Military Standards tables. 


(c) Quality-assurance process 

The quality-assurance process was introduced into manufacturing, align- 
ing production machinery with the requirements of product specifica- 
tions. ‘Control charts’ and, particularly in Japan, ‘quality circles’ were 
introduced. 
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(d) The quality-assurance concept 

Recognition of the possibilities opened up by quality assurance led to a 
reconsideration of the product in terms of feasibility of achievement; 
and this in turn led to decisions based on both quality and feasibility. 
The term ‘customer satisfaction’ came into use. 


(e) The ‘total quality’ concept 
Quality is now an important consideration both ‘upstream’ (marketing, 
production) and ‘downstream’ (sales, after-sales service): so the loop is 
closed. The concept first appeared in the USA with Feigenbaum (total 
quality control); this was an important breakthrough, all parts of the 
enterprise now being involved. 


(f) Achievement of total quality: ‘total quality-control system’ 
A total quality-control system involves bringing into play all the tech- 
niques that can affect the quality of the product. 

‘Meeting the customer’s needs’ means understanding his problems: 
there is often a world of difference between the stated need and the real 
need. It is marketing’s business to resolve this difficulty and draw up an 
appropriate specification. This is then studied by the research and 
development organization who decide on the technical methods to be 
employed and make a provisional forecast of the reliability of the 
product; further control of quality is exercised in the manufacturing 
process, and final control is at the stage of product release. The 
quality-assurance task is completed in the after-sales and maintenance 
services. 

It must never be forgotten that any enterprise is, above all, a 
collection of people: production will be affected by their health and 
quality of life; paying attention to these is a part of the quality process. 


1.3 WHAT IS TOTAL QUALITY? 


1.3.1 Definition 


‘Meeting the needs of, or providing the service required by, the 
customer or the user.’ 

Hidden in this statement are a number of points to which serious 
attention should be paid: 


@ the reliability of the product or service, 
@ the performance characteristics (of the product); 
@ its durability; 

@ its maintainability; 
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@ its security; 
© its effect on its environment (which must be acceptable); 
@ the ownership costs. 


These correspond to the AFNOR (the French official standards body) 
definition of quality; to extend the definition to total quality we must 
add consideration of the extent to which production of the product or 
provision of the service contributes to the satisfaction of the people 
involved in the enterprise: the shareholders and the staff. Thus, all parts 
of the enterprise are involved (Figs 1.3 and 1.4). 


1.3.2 Terms used in connection with quality 


(a) Quality assurance 

The accepted meaning of quality assurance is the laying down of a 
consistent set of standards and actions aimed at giving confidence in the 
achievement of quality. In practice this means compiling a quality 
manual and making sure that the engineer responsible for quality 
follows the key requirements in it. 


(b) Quality audit 
A quality audit is the detailed examination of 


e the product, 
@ the manufacturing process and 
@ the organization, in the context of quality. 


The existence of a quality standard is implied here. 


(c) Certification 

Certification involves the formal declaration by a recognized body that a 
product, a service or an enterprise meets a stated level of quality. It 
may take the form of the issue of a certificate or the authorization for 
the product to carry a particular label. 


(d) Marketing 
The role of marketing is crucial: for one thing it is upstream of 
manufacture, investigating where and in what volume the product will 
be sold, and for another it draws up the specification. A badly defined 
specification can involve the supplier in serious costs resulting from 
customer dissatisfaction and correction of errors. 

It is also the business of marketing to ensure that the customers are 
aware of, and appreciate, the quality of the product. 
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Study competitive Study customer 
d 
ogi i ie 
™—~ 
Draw up specification 
Prediction of EE 
reliabilit ‘~ 
y R&D ua Study customer 


reports 


Compare possible all 


technical solutions ™N Estimate life-cycle 


costs 
Manufacturing Control of 
Inward control Ly methods ™~ manufacturing process 
of raw materials aS 
Certification 


Study of 
customer's problems Nea, Service guaranteed 


After-sales service to customer 


Maintenance 


Creation of w Study of 
reliability database maintenance costs 


Figure 1.4 Activities required for achieving total quality. 


(e) Quality diagnosis 
Quality diagnosis is a matter of ensuring awareness of the key points 
concerning quality, which are as follows. 


Costs resulting from poor performance These are distributed among the 
various parts of the enterprise: 


@ internal — scrap, reworking, injuries, removal of pollutants etc. 

® external — costs of meeting guarantees and of providing after-sales 
services, penalties for late delivery etc. 

e detection of faults and out-of-course events; post-delivery checks and 
tests 


The overall cost of lack of quality is the sum of these separate costs. 
Technical backwardness This is the failure of the enterprise to keep 


abreast of technical developments — for example, numerical control, use 
of information technology. It can result in loss of market share. 
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Level of customer satisfaction A questionnaire can help to give an idea 
of the enterprise’s image in the marketplace, and of the return from the 
costs incurred in providing reliability, maintenance, performance etc. 


Overall quality can be expressed as the product of its components: 
O(overall) = Q(specification) x Q(design) x Q(execution) 
x Q(exploitation) 


(f) Ownership costs 

It is not sufficient, when considering the purchase of a machine or a 
product, simply to compare the cost-to-performance ratios of the various 
possibilities: maintenance costs must also be taken into account. This 
explains the importance of predictions of reliability as an aid to making 
the choice. 


(g) Quality circles, quality tools 
A quality circle in an enterprise is a voluntary group of six to eight staff, 
with a leader, who meet regularly to discuss and if possible to solve 
problems, technical and other, concerning quality. They will always have 
real-life situations in mind, while aiming at the ideal. 

Such a group will often be concerned with cause-and-effect investiga- 
tions and with priorities. 
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Reliability in the choice of 
technology 


The first essential is that the specification is fully defined, for on this the 
success of the project depends. This done, the reliability or ‘operational 
security’ of each possible technology must be investigated, so as to make 
possible the optimum choice. The term operational security (OS), which 
is coming more and more into use, means more than reliability and is 
understood as including 


e reliability — ability to work without failure 

@ maintainability — ability to be restored quickly to working condition 
after failure 

®@ availability — being in working condition when required 

® security — remaining safe in case of failure 


2.1 QUANTITATIVE ANALYSIS: FAILURE RATES, 
RELIABILITY LAWS 


The instantaneous failure rate A(t) is defined by the statement that the 
probability that the device under consideration will fail in the (infinite- 
simal) interval (t,t + dt), having operated without failure up to time f, 
is A(t) dt. The cumulative failure function F(t) is the probability that 
the device has failed at least once before time t is reached, and the 
reliability function R(t) is the probability that it has not failed up to this 
time, i.e. that it has operated reliably. Clearly 


R(t) =1- F(t) 


F(t + dt) is the probability that the device fails at least once up to time 
t+ dt, and this will happen either if it fails not later than time ¢ or if it 
does not fail up to that time and fails in the interval (t,t + dt). This 
gives the relation 


F(t + dt) = F(t) + R(t) A(t) dt 


Quantitative analysis: failure rates, reliability laws 11 
i.e. since R(t) =1— F(4),; 
Et seadt jean h UL) 
de) 
aa SORE) 
~ 1- F(t) 


If we reckon time from t = 0, F(t) = 0 at t = 0, and so integrating from 
0 to t we have 


A(t)dt = 


[a ae = Tih oEeohe a) 
and hence 


R(t) = exp - [aw ar 


F(t) =1- exp|— [400 ar 


A third function, the failure probability density function f(t), is 
dF 
f(t) = ane ure exp |- [, Ar) ar 
cea 


The equations between R(t), F(t), f(t) and A(t) are the most general 
expressions for the laws of reliability. 

An important quantity is the average time of fault-free operation, or 
the mean time between failures, MTBF. This is the mathematical 
expectation of the time to fail: 


MTBF = f t f(t) dt 
If we integrate this by parts we get the equivalent expression 


MTBF = I, R(t) dt 


2.1.1 Reliability models 


(a) Constant failure rate: the exponential law 

In general, electronic components that have reached a state of maturity 
show a constant failure rate; this is expressed by putting A(t) equal to a 
constant, say A, so that 
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[a dr=At R(t) = exp(— At) 


Example 1 

If A=2%x10~° failures per hour and t=500h, then At =0.001. So 
R(t = 500) = exp(— 0.001) = 0.999 and MTBF= f~exp(-At)dr= 
1/1 5 x 10° h. 


(b) The log-normal model 
The log-normal model gives a good representation of mechanical fatigue 
or wear. The probability density function is 


1 1/Int — m\?] 1 
I) aaa exp|- 4 0 kh 


where m and o are the mean and standard deviation respectively of Inf. 

Calculations with this model are carried out most easily in terms of 
the reduced variable u = (Int — m)/o, which is distributed normally with 
mean zero and standard deviation unity. Making this substitution we 
find 


MTBF = |, rf(¢) dt = exp(m + 40’) 


Example 2 

The lifetime of the con-rods of an automobile engine follows a log- 
normal law with parameters m=5, o =1.4, time being measured in 
hours. Find (1) the reliability after 300 h and (2) the MTBF. 


(1) wu = (1In300 — 5)/1.4 = 0.502. From Table1 (p.154) we find that 
F(u) = 0.692, and so R(t = 300) = 1 — F(0.502) = 0.308 (which is 
poor). 

(2) m+ 50° = 5.98, and so MTBF = exp (5.98) = 395 h. 


(c) A more general law: the Weibull model 

(i) General form of Weibull’s law 
The most general form of Weibull’s law includes many of the simpler 
models and is given by 


RG@)osd eves exp|-(“—2) 


The instantaneous failure rate A(t) is 
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Rw 


Figure 2.1 Examples of Weibull law. 


yea"! 
RGA Ain 
B is the shape parameter, n is the scale parameter and y is the location 
parameter. A study of the form of A(t) shows that 


for B<1  A(Z£) is a decreasing function of t 
for 8 =1 A(t) is constant, equal to 1/n 
for 6B >41 A(t) is an increasing function of t 


For the particular case B=1, y=1, Weibull’s law reduces to the 
exponential law: 


R(t) = exp(— t/n) 


i.e. the exponential law with parameter A = 1/n. 
For B 23 Weibull’s law approximates to the normal law, more closely 
as f increases. 


(ii) Estimation of the parameters 
A basic problem in connection with the use of Weibull’s law is 
estimating the values of the parameters 6, n and y for given data. Two 
methods are available: 


@ a purely numerical method, which leads to differential equations that 
are difficult to solve and consequently is little used; 

© a graphical method which uses special paper, called Weibull paper, 
ruled with functional scales. This is the method most used, and the 
one we now describe. 


The scales are as follows: ordinate, Y = InIn {1/[1 — F(t)]}; abscissa, 
X =Int. 
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The case y = 0 corresponds to the assumption that the origin of time 
is known and is given by the data. Then 


7 Ve 
R(t) =1- Fi) = exp|-(“) | 


: =F a5) 


[- Bint — Blnyn 


and so 


Y= taal — FO) 


Xo— nee 


X and Y are both functions of t; but if we put A = 6, B = 6 lnn, which 
are constants, we have 


Y= AX -B 


which is the equation of a straight line. Thus if a set of observations of 
F(t) at a number of values of ¢ are described by a Weibull law with 
y = 0, then these should lie on a straight line on Weibull paper. 

Figure 2.2 illustrates this. On this paper the origin for Y is the 
ordinate F(t)=0.632 (or 63.2%), because if Y=0 then 
In {1/[1 — F(t)]} =1 and so F(t) = 1 —- 1/e = 0.632.. 

The shape parameter # is the slope of the line. To find its value we 
draw a line through the point (t=1, F(t) =0.632) (marked on the 
paper) parallel to the line on which the data points lie and read the 
value at the intersection of this with the # scale. 

The value of the scale parameter 7 is read at the intersection of the 
data line with the parallel to the X axis through the ordinate 
F(t) = 0.632, because there Y = 0 and so AX — B =0, which from the 
definitions of A and B gives X =Inyn. Since X =Int, this point is 
mn 

In the example of Fig. 2.2, 6B = 1.5 and n = 20000 h. 

For the case y = 0 which we are considering 


MTBF = Me R(t) dt = nT + 1/8) 


where I" denotes the gamma function: formulae and a table are given in 
Appendix 9. For the above example we find, for f=1.5, 
(1 + 1/8) = 0.9027, and so (to a realistic accuracy) 


MTBF = 20000 x 0.9027 = 18000 h 


For y>0 the data cannot be linearized by the above process; instead, 
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the points will lie on a curve that has a vertical asymptote, and y is 
given by the value of t at which this asymptote intersects the t axis. This 
follows from the fact that, when t=y, F(t)=0, and so 
Y =InIn1 = 1n0 = —o. Figure 2.3 illustrates this. 

To find the parameters in this case we first estimate the position of 
the asymptote, as in Fig. 2.3, and obtain a first estimate for y, say y’. 
We now change the ¢ scale to t’ =t-— yy’ and repeat the linearization 
process with the new time scale. If this gives an acceptable approxima- 
tion to a straight line then y’ is a sufficiently good approximation to the 
true value of y and we can continue as before, finding f and 7. If not, 
we go through the estimation process again, estimate a value y” and plot 
the data again with the time scale t” = t' — y’ =t — y’ — y”; and so on. 
But if a third repeat of this process does not give an adequately linear 
plot we must conclude that the data do not follow a Weibull law; they 
may follow a mixture of Weibull laws with different parameters, or 
some quite different law. 

For y <0 we have t— y>0 for t>0 and F(t) 1 — exp[-(— y/n)] 
as t—0. Thus y (= InIn1/{1—- F(t)])—fln(—y/n) as X 
(= Int)— —, i.e. the curve of Y against X has a horizontal asymp- 
tote. 

One way of proceeding in this case is to try a succession of estimates 
for y until an acceptably linear plot is obtained and then to continue as 
before. Figure 2.4 illustrates this. 


99.9% 


63.2 


0.1 


Estimation of y (for y < 0) 
Figure 2.4 Fitting Weibull law (y <0). 
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(iii) Another method for determining y 
Using the same change of variables as before, an estimate for the value 
of y can be computed as follows: 


Vg ma Pi (CG ca Gy) 
CX vat rT X ju)) a kX on ai ae 
where X nx is the value of X corresponding to the maximum value Y jax 
of Y, Xmin is the value of X corresponding to the minimum value Y jin 


of Y and X,, corresponds to the midpoint Y,, between Y,,,, and Y min 
measured on a linear scale. This is illustrated in Fig. 2.5. 


y= Xn 


0 Ximin. Am Xmax. 


Figure 2.5 Weibull law: Estimation of y using method in 2.1.1 (cxiii) 


Example 3 
Given that the following data can be represented by a general Weibull 
law, find the values of the three parameters Y, 6 and 7. 


lifetimes in hours: 
705, 812, 902, 995, 1070, 1171, 1301, 1440, 1650 


From these values 


Ymax = 0.9° > Xn =1650 
Yin = 0.1 S X..= 705 
Yn =038 +X, = 980 
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Substituting in the above formula we obtain 
y = 980 — (1650 — 980)(980 — 705)/[(1650 — 980) — (980 — 705)] 


= 513h 


Knowing y we can linearize the curve and find the other parameters; we 
obtain the values 6 = 1.8, 7 = 700 h. 


Example 4 
The lifetime values (h) of a mechanical system are as follows: 


5, 112, 202, 295, 370, 471, 601, 740, 905 


The corresponding pairs of values of Y and X are now (0.9, 905), 
(0.1, 5) and (0.38, 275), and the same calculation as in Example 3 gives 


y=197h BH=18 =705h 


(iv) Estimation of F(t) 
There are two methods for estimating F(t;). 


(1) Method of median ranks (for small samples): 


F Ws 3 nN; as 0.3 
ks n + 0.4 
There are tables of this (see pp. 170-177). 


(2) Method of mean ranks (more commonly used): 


(v) Application: an example 
A sample of nine ball bearings has been put into service as a test of a 
new production series. The results were the following lifetimes (h): 


801, 312, 402, 205, 671, 1150, 940, 495, 570 


(1) Assuming a Weibull law, find the parameters. 

(2) Compute the MTBF. 

(3) Find, both graphically and by computation, the reliability after 
600 h. 


(1) We start by putting the observations in increasing order and 
tabulating the distribution function 
i i 


F(i) = 


(Table 2.1). We next plot these on Weibull paper with lifetimes as 
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abscissae and F(i) as ordinates; this gives the graph in Fig. 2.6. From 
the figure we find y=0 (because the points lie on a straight line), 
B=1.8 and n= 710h. 

(2) MTBF = E(t) = 7I(1 + 1/8) = 710 x T'(1.555) = 631 h from Table 
Sick: 

(3) Computationally, R(t = 600) = exp[—(600/710)'*] = 0.480. From 
the graph, at t = 600, F(t) = 0.52, and so R(t) = 1 — F(t) = 0.48. 


Table 2.1 Observed data for Figure 2.6 


Position Lifetime F(i)% Position Lifetime F(i)% 
1 205 10 6 671 60 
2 312 20 y | 801 70 
3 402 30 8 940 80 
4 495 40 9 1150 90 
5 570 50 


99.9% 


63.2 


F(t) 


0.1 


Figure 2.6 Weibull law: Fitting data from Table 2.1. 


(d) Mixture of Weibull laws (use of method of median ranks) 
A mixture of different Weibull laws may be needed to model a set of 
failure data because 


© the items may have come from different populations 
@ several different modes of failure may coexist simultaneously 
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This is illustrated in Fig. 2.7 which shows three different populations 
coexisting. 


100% 


63.2 


Figure 2.7 Weibull Law: Three populations coexisting. 


(i) Solution of the mixed law problem 
The data (times of fault-free operation) are plotted after applying the 
method of median ranks, 


2) n; — 0.3 

n+ 0.4 
(or using Johnson’s table, for small samples). 

The size of each subpopulation is found by counting the number of 
items that belong to it, say n,, nz, n3, and an estimate of the proportion 
of each in the total population is given by P,=n,/n etc. Each 
subpopulation is studied independently, either by plotting on Weibull 
paper or by computation; thus for the first we would have 
> ny — 0.3 


lA 


F(t,;) = 


and similarly for the others. 
The reliability model for this system with three subpopulations is then 
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The following lifetimes are observed, given in ascending order: 235, 390, 
540, 580, 730, 766, 800, 850, 900, 940, 980, 1100, 1150, 1200, 1240, 
1310, 1400, 1455. Find the model. 


Ri23(t) = SS P;exp 


Plotting on Weibull paper indicates that there are two populations P, 
and P,. For P,; we have the values given in Table 2.2 and from the plot 
y =0, B=2.3 and n = 450. 


Table 2.2 Values for population P, 


Order Lifetime F(t;) 


1 PB) 20.5 
2 390 50 
3 540 79.41 


For P, the values are given in Table 2.3, and y=0, B=4 and 
7 = 1020. Hence 


3 ¥\23)" 46 fe 6 
R(t) = me|- (=) ry ge|- aan) | 


Table 2.3 Values for population P, 


Order Lifetime F(t;) Order Lifetime F(t;) 
1 690 5.88 9 1020 52.94 
2 730 11.76 10 1100 58.82 
3 766 17.64 ity 1150 64.70 
4 800 23,53 12 1200 70.58 
5 850 29.41 13 1240 76.47 
6 900 35.29 14 1310 82.35 
7 940 41.17 15 1400 88.23 
8 980 47.06 16 1455 94.12 


(e) Program for finding the Weibull parameters, assuming that y = 0 
The program uses the transformations 
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1 
X, = Int, = inn 
; = Int, va nin| | 
Estimates of the parameters B and n (actually, n®) are found by 
least-squares fitting of a straight line (Fig. 2.8). 


10 CLS 

20 DIM T(100) 

30 PRINT “WEIBULL MODEL” 

40 PRINT “ENTER TBF ONE AT A TIME” 
50 PRINT “DO TBF=0 TO PROCESS DATA” 
60 I=0 

70 I=I+1 

80 PRINT I 

90 INPUT “ TBF=”;T(I) 

100 IF T(I)=0 THEN 120 

110 GOTO 70 

120 N=I 

130 I=0 

140° T=0 

150 F 
160. T2= 
170 F2 
180° F5= 
190 PRINT “number of TBF taken=” 
200 N=N-1 

210 PRINT N 

220 N=N+1 

230 I=0 

240 FOR I=1 TO N-1 

250 T=T+LOG(T(I)) 

260 F1l=1/(1-I/N) 

270 F3=LOG(F1) 

280 F=F+LOG(F3) 

290 T2=T2+LOG(T(I)) 2 

300 F2=F2+LOG(F3) 2 

310 F4=LOG(T(I)) «LOG(F3) 

320 FS=F5+F4 

330 NEXT I 

340 N=N-1 

350 B=(F5—T*F/N)/(T2-T 2/N) 
360 E=T/N-F/N/B 

370 E1=EXP(E) 

380 R1=F5—T+F/N 

390 R2=(T2—T 2/N)+*(F2—-F 2/N) 


Figure 2.8 BASIC program for computing Weibull parameters (y = 0). 
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400 R3=SQR(R2) 

410 R=R1/R3 

420 CLS:PRINT “RESULT OF PROCESSING:” 

430 PRINT“BETA=”;B 

440 PRINT“ETA=”;E1 

450 PRINT “LINEAR CORRELATION :R=”;R 

460 PRINT ‘‘Compute MTBF if YES enter 0 if NO 20” 
470 INPUTG 

480 IF G=0 THEN GOTO 490 ELSE GOTO 570 

490 X=0 

500 P=.001 

510 X=X+P 

520 G=G+ABS(LOG(1/X)) (1/B) *P 

530 IF X>=1 THEN 550 

540 GOTO 510 

550 X=E1*G 

560 PRINT“MTBF=”;X 

570 PRINT ‘Compute the R(T) and F(T)” 

580 INPUT “‘to compute R(T) and F(T) enter t if NO do 0”; T3 
590 IF T3=0 THEN GOTO 670 

600 R1=(T3/E1) B 

610 R2=EXP(—RI1) 

620 PRINT “For t=”;T3;“R(T) is equal to’’;R2 

630 R4=1-—R2 

640 PRINT “For t=’’;T3;‘‘F(T) is equal to’’;R4 

650 INPUT “To repeat calculation enter t if NO 0”’;T3 
660 IF T3=0 THEN GOTO 670 ELSE GOTO 590 
670 END 


Figure 2.8 (continued) 


2.1.2 Verification of the models 


Any model constructed in a reliability study will be based on a sample 
drawn from the population being investigated, and some assumption will 
always be made concerning the distribution law for that population — 
exponential, log normal etc. There is therefore the question of the 
validity of this assumption, and this can be answered by applying what 
are called goodness-of-fit tests. In using these statistical tests we must 
always recognize that there is a risk of being wrong, measured by the 
probability a that the test will give the wrong result. @ is called the 
significance level of the test, and we aim to make its value small. 


(a) The chi-squared (32) test 
The condition for the x test to be applicable is that there are at least 50 
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observations: n > 50. It is usual to group the observations into classes so 
that there are at least five in each class; the classes need not be at 
regular intervals. 

The test is based on the differences between the numbers of observa- 
tions in each class and the number predicted by the model; the measure 
of this difference used by the test is 


where r is the number of classes, n; is the number of observations in 
class i, n is the total number in the sample (= 2 n,), p; is the probability 
that an observation will be in class i and np; is the expected (theoretical) 
number in class i. E is distributed approximately according to the y* law 
with v degrees of freedom where v= r — k — 1 and k is the number of 
parameters whose values have to be estimated in deriving the model. 
This depends on the underlying law assumed: for example, k = 1 for the 
exponential law, k = 2 for the normal (Gaussian) law and k = 3 for the 
Weibull law. 

y° is a function of two variables, the degrees of freedom v and the 
significance level aw; for given w 


Bl Ne Soiedbaaa’ aah 


and the test is that if E > 7, _, we reject the hypothesis on which our 
theoretical model is based. 

It should be noted that other tests can be used, e.g. the Kolmo- 
gorov-Smirnov test (see Lyonnet, Maintenance Planning, Chapman & 
Hall). 


2.1.3 Predicting reliability 


In order to estimate the cost of a maintenance service and to decide 
how to implement this, it is important to be able to predict the 
reliability of the equipment that is to be maintained; this is also 
important for the choice of techniques to be used. Published data are 
available on which the prediction can be based. 


(a) Relevant databases 

There are tables relating to currently used electronic and mechanical 
components; those most commonly used are the following: in France, 
those published by CNET (NPRD 1 and 2) and by EDF; in the USA, 
Rome Air Development Center (RADC), NASA, US Navy 
(FARADA) and AVCO Corporation. The tables give 
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the name (identifier) of the item, 

the MTBF, 

the failure rate, either average or calculated on the assumption of a 
constant A(t), 

basic statistical information (e.g. confidence intervals), 

a multiplier to be applied to the given failure rate when the 
equipment is used under each of a number of stated conditions. 


Two comments are relevant. 


. There are fundamental differences between electronic and mechanical 


systems. 

For electronic components the statistical information is more 
important than for mechanical components and the MTBF derived 
from this is more reliable. The failure rates are usually constant and 
can be taken as the values given in the tables. 

This does not hold for mechanical components. In practice, failure 
rates are found not to be constant, mechanical components are less 
well differentiated than electronic components and there is less 
statistical information. 


. The conditions under which the equipment is actually used are often 


very different from those assumed in the tables. The multipliers that 
should be applied in the various environments do not always take 
into account 


the installation conditions 
vibration 

temperature 

dust 

corrosion 

mechanical constraints 


The general conclusion is that results obtained on the basis of the 
published tables should be treated with caution, especially in the case of 
mechanical systems. 


For illustration, Table 2.4 gives an extract from the CNET NPRD-2 


and RADC-NPRD-3 tables. 


(i) The RADC NPRD-3 tables 


These provide a databank for mechanical systems and reliability figures 
are included. The headings have the following meanings: 


CLASS a family of components having the same 


function 


TYPE the particular member of that family 
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ENVIRONMENT coded as follows 


DOR Dormant: the item is connected to the system 
but is out of operation for long periods 

SAT Satellite: in orbit around the earth; no access 
for maintenance 

GRF Fixed terrestrial installation; permanent; 
ventilated; maintenance by military personnel 

GRM Mobile terrestrial installation; conditions 


harsher than for GRS — vibration, shocks; 
maintenance more difficult 


(b) System reliability 
From the point of view of reliability the aims in constructing a system 
made up of a number of components are as follows: 


® to satisfy the customer’s requirements as expressed in the reliability 
specification, or his needs if not so specified; 

®@ to choose an appropriate technology, using the reliability cost ratio as 
criterion; 

® to improve the reliability by bringing to light the critical points; it 
should be possible to make a prediction of the reliability in the design 
stage. 


A system consists of a set of elements or subsystems each of which 
provides one or more stated functions; thus the design proceeds by 
breaking the system down into elements, for each of which a numerical 
value for the reliability can be given, and then constructing a representa- 
tion of the organization of these together to form the complete system. 
This is called constructing a block diagram for the system. 


(i) Block-schematic reliability calculation 
A series system (Fig. 2.9) fails if any one of its components or subsys- 
tems fails; if R,(t) is the reliability function for component or subsystem 
i, and if all are independent, the system reliability R,(t) is 


ROS Rit) Ra) & >. & RAG) STR) 


TH} 1}---a-— 


Figure 2.9 Block schematic for serial system. 
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A parallel system (Fig.2.10) fails if and only if every one of its 


components or subsystems fails; since failure corresponds to the failure 
function F(t) = 1, the equation now is 


Fp) = ILE) 
or, since R(t) =1— F(t), 


R,(t) = 1—- fl — R(t] 


Figure 2.10 Block schematic for parallel system. 


A general system can always be represented as a collection of series 


and parallel subsystems, themselves connected in series and/or parallel. 
Thus for the system in Fig. 2.11 


Figure 2.11 Block schematic for series parallel system. 


Ry = R,R,[1 his R;)]{1 athe R,)] 
Kae Ce a) (ke 
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(ii) Reliability scenario 
An item of equipment may operate in different modes, in different 
places and under different conditions at different times — for example, 
idle, transported to a new site, put into operation there. The calculation 
of its predicted reliability must take such a possibility into account. Thus 
if the item under consideration, at the end of a period ¢, will have spent 
times t,, ft and ft; in modes 1, 2 and 3 respectively, the reliability is 


R(t) = R(t1)R(t2)R(ts) 


Example 6 
A piece of radar equipment is placed for 2 h per day on the bridge of a 
ship and for 22 h on shore at sea level; if the failure rates are 


on shore, sea levelA, = 2 X 10~° failures per hour 


on board A, = 3.6 X 10~° failures per hour 


what is the reliability after 500 days? 


Here t,(on shore) = 11000 h and t¢,(on board) = 1000 h; therefore 
Ait, = 0.022 At, = 0.0036 
R(t = 500 days) = exp (—0.022) exp (—0.0036) = 0.9747. 


Example 7 
A machining unit has four machines organized as in Fig. 2.12. If the 
separate reliabilities are 


Ra — 0.95 Rp = Rp — 0.97 Rec = 0.98 


the overall reliability is 


R, = (0.95)[1 — (1 — 0.97)3](0.98) = 0.93 


Figure 2.12 Block schematic for reliability calculation. 
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(iii) Calculation of the MTBF (or MTTF) 
MTBF (mean time between failures) is relevant for repairable systems, 
MTTF (mean time to failure) for non-repairable systems. 

The general result (see p. 11) is 


MTBF = I, R(t) dt 


Evaluation of the integral is particularly simple when the failure rates A; 
are constant, i.e. when the exponential law R(t) = exp (—At) applies. 

For a serial system, integrating the expression R,(t) = I1,R,(t) gives 

1 1 
MTBF = SF ay 

if all the A; are equal. 

For a parallel system, integrating R,(t) = 1 — II,[1 — R,(t)] gives 

1 1 1 

MTBF =) 3 SS De 
Xz rey ink Aiea eh iAg 


i 
IF] l#jJ#k 


If all the A; are equal the expression is 


ff ¢ 1 182 1 orl 
MTBF = [, (1 - [1 - exp(-any*jar=F(14545+...+ 4} 


[Translator’s comment: this result shows first that, provided that all the 
elements are truly independent (very important), the MTBF can be made 
as great as one wishes, because the series 1+ 1/2+1/3+... diverges; 
but second that there is a law of diminishing returns, and the gain from 
adding a further element in parallel decreases steadily. Thus the MTBF 
can be doubled by putting four elements in parallel, but to multiply it by 
3 needs 11.] : 


(iv) Interrupted tests 
There are two main types: 


Sequential in which the test is stopped after some agreed number of 
faults, say C, have been recorded 
Truncated in which the test is stopped after an agreed time, T say 


Both types can be conducted with or without replacement. 

Suppose that a sequential test without replacement starts at t = 0 with 
n items working and that the successive failures occur at instants 
t,, t>,..., t, when c failures have been recorded and the test is stopped. 
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The total amount of fault-free working time is then 
T= Trak oe ee eee 

and the estimate of MTBF is 


fie. Sh = eril) 


(6 


MTBF = 


In sequential testing with replacement T = nt, and MTBF = nt,/c. 

The truncated test without replacement finishes at a time ft; agreed in 
advance; suppose k failures are recorded, occurring at instants f,, 
to, .. +5 ty (te < ty). The total good time is 


ISR aoe on hg 


and MTBF = T/k. 

For the truncated test with replacement, using the same notation, 
T = nt, and MTBF = T/k. 

As we shall show in Chapter 6, we can construct a confidence interval 
for each of these estimates. Further, the MTBF can be used to estimate 
the corresponding value of A. 


2.1.4 Markov chains: reliability and availability 


Markov chains are a mathematical technique that enables us to compute 
the reliability of a system. A system consists of a set of elements 
connected in series and/or parallel, as in Fig. 2.13. At any instant it will 
be in one of a number of possible states and may or may not change to 
another state. We make the following basic assumptions. 


Figure 2.13 Series parallel system. 
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1. The possible states are numbered in such a way that when the system 
is in state 7 it can change only into i — 1 or i + 1; this means that the 
state it is in at any instant depends only on the two neighbouring 
states. 

2. With the state changes corresponding to the failure of an element or 
the repair of a failed element, both failure and repair times follow an 
exponential law with constant rates A and wu respectively. 


(a) Transition graph and equations 

We consider a system of n elements. We define the state i as that in 
which i elements are working satisfactorily; thus in state n the whole 
system is fault free and in state 0 it has failed completely - a 
breakdown. We denote the probability of changing from state i to state 
] by pj; it follows from (1) above that j can be only i— 1 ori +1. 

It is convenient (and illuminating) to represent the state changes by a 
labelled directed graph; Fig. 2.14 illustrates this for a three-element 
system. Let P(i, t) be the probability that the system is in state i at time 
t. It will be in state 7 at time ¢ + dt if 


@ it was in state i at time ¢ and did not change during the interval drt or 
® it was in state 7 — 1 at time ¢ and changed to / during dt or: 
@ it was in state i + 1 at time ¢ and changed to i during dt. 


© e 


p01 p23 


Figure 2.14 Transition graph. 


These are the only possibilities; since they are independent we can add 
the probabilities: 


PG it + dt)-= PG, 1) prygdd) = piajdt) +oP Ga 1pt) p;-r7dt 
POH, mde 


Simplifying, we obtain © 
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Plisterdt)cawP@) 
dt 


= PU.+ 1, Ping + PU 1 picu 
— P(i, t)(pis1i + Pi-1,,) + terms of order dt 


Letting dt — 0 
dP(i, t) 


ain PWT, Paty ee — 1p Ete Opa: 1 eee 


This is the most general state-change equation; assuming that the system 
is in full working condition at t = 0 the initial conditions are 


Pi, 0) =Ofori#n P(n, 0) = 1 


There is also the condition that at any time ¢ the system must be in one 


or other of the possible states 0, 1, 2, ..., n; thus 
DM Ge 
0 
Lee 
SP, 7PrA0 
0 


where P’ = dP/dt. The equations are linear differential equations for 
the probabilities P(i, t), and under the conditions that we have assumed 
the coefficients p; (the transition probabilities) are constant; they can 
therefore be solved by the Laplace transform method, as we shall show 
below. 

The equations can be written very easily with the help of loops added 
to the transition graph: at each state 7 we add a loop in which we write 
the sum of the probabilities p,; of changing to another state j with the 
signs reversed. Then dP(i, t)/dt is equal to the sum of the products of 
the transition probabilities associated with each arc arriving at state i by 
the state from which that arc started. This is illustrated in Fig. 2.15; 
applying the rule gives the equation for dP(i, t)/dt obtained above. 

To apply these equations to any actual problem we have to know the 
transition probabilities p,. For the reliability problem we are studying 
we have 


failure rate A = 


repair rate u = 


MTBF MTTR 


If there are a number of repair stations working independently the mean 
time to repair is reduced by that number and therefore w is increased. 
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p(i+1, 1) 


-p(i-1, 1) 


p(i, i+1) p(i-1, i) 


Figure 2.15 Transition graph leading to state-change differential equations. 


(b) Availability of a simple system: solution of the equations 
The simplest system consists of a single element; this element could of 
course be the representation of a more complex system for which we 
knew the failure and repair rates for the system as a whole. 

There are now two and only two states, i=0 or i=1, with the 
Markov chain as in Fig. 2.16. The equations are 


P'(1, t) = —AP(, t) + uP’(0, 2) 
P’'(0, t)= AP(, t) — uP, t) 
with 
P,0)=1 #£P(O,0)=0 THE LIBRARY | 
FORTH VALLEY COLLEGE 
WESTSURN AVENUE 
nN FALKIRK FiK4 SST 


STIRLINGSHIRE 
Welephone (0324) 63500 


Figure 2.16 Markov chain for a two state system. 


Applying the Laplace transform to these equations (L; is the trans- 
form of P(i, t)) 


pl, a Pd, 0) = —AL; ae Lo 
pLo = P(O, 0) = AL, = Lo 
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Putting P(1,0) =1, P(0,0) = 0, we have in matrix form 


pak = ie ey | 
=—4 p+ p} \Lo 0 
We are interested only in the availability of the system, i.e. in the state 
in which everything is working, P(1,t). Therefore we need only 
consider the solution for L;: 
= u+ p 
put A+ p) 


To recover P(1, t) by the inverse transform L — 1 we must put this into 
partial fractions 


Ly 


A B 
hy See 
Dp fico eee p 
We find 
A 
Ae ee 
(abi aae ae 
We now have 
A B 
Le = + 
Dp PEED 


where a = w+ A, and so from tables of the Laplace transform we find 
P(1, t) = A + Bexp(—at) 


P(1, t), the probability that the system is in a working state, is the 
availability A(t); thus we have 


u A 
A(t) = . = 
(t) ae oe ee (u + A)t] 
As t increases, A(t) tends to the constant value u/(u+ A). Thus the 
result is that the availability of the system approaches the steady state 
value 


u MTBF 


u+A MTBF +MTTIR 


(c) System with redundant elements 

Suppose there are n identical elements, each with constant failure and 
repair rates A and u, connected in parallel; the transition graph is given 
in Fig. 2.17. Using the shortened notation P; for P(i, t), the equations 
are 
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PSAP 
Pi =(i+-1)AP;,,-—GAt+ yp) P;+uP;-1 i=h2,...n-1 
Po = AP, — wPo 


with initial conditions P,=1, P;=0, i#n at t=0. These can be 
solved by the same method as before, but unless n is very small it is 
advisable to use a software package. We consider here the case n = 2, 


Pi = =2AP st uP, 
Pi = JAP, — (ue ANP, 
Pos AP, — UP 


with P,; = 1 and P,; = Pp) =0 at t=0. 


Figure 2.17 Markov chain for multi-state system. 


Transforming and collecting terms as before, we obtain the equations 
in matrix form: 


PPA = 0 B 1 
=) pene =i LE, 7=T70 


Solving, inverting and taking the constant term we find 
w+ 2ua 


(ee ee Eas 
re calinlky TiiRa Bie 
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The graph of A(t) as a function of ¢ is shown in Fig. 2.18. 


yw? + 2074+ 2ur 


0) 


Figure 2.18 Availability function for Markov system. 


Results such as these enable rational decisions to be taken on such 
things as the choice of technology, the amount of redundancy to build 
into ihe system, the number of repair stations, the amount of replace- 
ment stock to carry etc. It is advisable to use computer aids when 
dealing with complex systems. 


2.1.5 Simulation: the Monte Carlo method 


A simulation of a system enables a range of possibilities to be studied 
and hence an optimum situation to be defined. The so-called Monte 
Carlo method can be applied in this way to study reliability. 


(a) Principle of the method 

If the cumulative failure distribution function for the system is F(t) the 
method is based on the idea of choosing a random sample x from a 
population distributed according to F(x). In the reliability study x is a 
value for the lifetime of the element or system, say ¢; for the ith sample. 
Suppose we draw N samples and that N, is the number of these with 
t; > t;, where f, is the required time of fault-free operation. Then 


R(t.) = N,/N 


is an estimate of the reliability at f,. 


(b) Procedure 
There are five stages. 
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1. Obtain the distribution function f(t) for the lifetime of the equip- 
ment under investigation. 

2. Derive from this (by integration) the cumulative failure distribution 
function F(t) (the probability that there will be at least one failure by 
time f). 

3. Obtain a set of random numbers uniformly distributed between zero 
and unity. 

4. Construct a random sample of lifetimes as follows: choose at random 
one of the numbers in (3), r say, and from the graph or table of F(t) 
find the value of ¢, t, say, such that F(t,) =r. Repeat this until a 
sample of the required size, N say, has been drawn. 

5. Use the sample to estimate the reliability, as above. 


Uniform law 


Figure 2.19 Drawing random numbers distributed between 0 and 1. 


(c) Some applications 
Example 8 
The lifetime distribution for a gyroscopic system is 


t 
My) = ra exp(- as 
(t in hours). Find the reliability at t = 2000 h. 
F(t) = [ fii)dt =1- exp - al 


Ten numbers r are drawn at random from a set uniformly distributed 
between zero and unity, and Table 2.5 is constructed. Three of the ten 
trials give ¢ > 2000 and so the reliability is R(t = 2000) = 0.3. 
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Table 2.5 Monte Carlo solution for Example 8 


Drawing no. r= F(t,) by t, > 2000 
1 0.43793 927 No 
2 0.07496 124 No 
3 0.17405 306 No 
4 0.80966 2654 Yes 
5 0.65989 1725 No 
6 0.55400 1292 No 
7 0.72301 2054 Yes 
8 0.36504 a2. No 
9 0.00187 2 No 

10 0.90375 3745 Yes 


For comparison, calculation by the exponential law gives 0.286. 

There are in fact two random variables involved in a problem of this 
type: the lifetime of the equipment (t, above) and the time f¢, actually in 
service. If we draw a second set of samples to simulate the time in 
service we can compare the results with those for the lifetime to get 
another estimate of the reliability. If for convenience we rename tf, as f;, 
the new estimate is 
ndmiber (77 54,) 

Example 9 

For the equipment of Example 8, suppose the time in service is 
distributed normally with mean 2000 h and standard deviation 150. Find 
the new value of the reliability. 


The sampling for in-service time is done in the same way as for the 
lifetime, with the difference that F(t) is now the cumulative normal 
(Gaussian) distribution with m = 2000, o = 150; we obtain the required 
values from the table in terms of the reduced variable u = (t — m)/o, 
and having found a value u we convert this back to ¢ using t = ou + m. 

Thus if we choose the random number 0.9408 we find u = 1.56, i.e. 
F(1.56) = 0.9408, and so the value of ¢t is 150 x 1.56 + 2000 = 2234 h. 
The two drawings of random numbers (RNs) are given in Table 2.6, 
giving again the estimate R(t) = 3/10 = 0.3. 

This is a better procedure because it uses more values for the sample. 
Random numbers are easily generated by computer, and this has led to 
the development of software for simulation. 
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Table 2.6 Monte Carlo solution for Example 9 


RN(1) ty RN(s) fs Spits) 
1 0.43973 927 0.94080 2234 No 
2 0.07498 124 0.27777 1911 No 
3 0.17405 308 0.09621 1804 No 
4 0.80968 2654 0.45577 1982 Yes 
5 0.65989 1726 0.78282 ale lgh No 
6 0.55400 1292 0.10039 1808 No 
7 0.72301 2054 0.19572 1872 Yes 
8 0.36504 727 0.09306 1802 No 
9 0.00187 3 0.89518 2188 No 
10 0.90375 3745 0.900041 2193 Yes 


2.2 QUALITATIVE ANALYSIS 


2.2.1 Use of failure mode analysis for quality improvement 


Failure mode analysis (FMA) is a rigorous procedure for detecting 
potential faults, in which both the probability of the fault’s occurring 
and the seriousness of the situation should it occur are taken into 
account. It is a very valuable tool for reducing the risk of equipment 
operating badly or failing in service and should be included in any policy 
for total quality control. 

FMA can be applied to a product or to a process. In a manufacturing 
industry such as the automobile or the aeronautics industry, that uses a 
number of subcontractors, the main contractor should require all the 
suppliers and subcontractors to use FMA for products and procedures in 
order to guarantee quality. 


For products: all possible failure modes of the system or 
subsystem that is being designed are noted and are 
taken into account in the analysis. 

For procedures: the analysis takes account of all failures that can 
result from the manufacturing processess — 
assembly, casting etc.; the research and 
development organizations are all involved here. 


2.2.2 The practice of ‘Product FMA’ 


In implementing a policy of total quality control the analysis must be 
applied to 
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@ all new components 
e all components that have been modified or are to be used in new 


circumstances 


(a) Failure analysis 

Failure analysis is the most difficult part, and the one that demands the 
most skill. It involves extending the design calculations done by the 
research and development department to take account of all the 
influences that could bear on the components. Thus failures could be 
caused by 


deformation 
fatigue 

crack propagation 
brittle fracture 
vibration 
seizing-up 
leakages 
corrosion 
short-circuits 


(b) The aim of FMA: customer relations 

The prime importance of FMA is that it brings to light the critical 
issues, so that either the possibility of critical situations arising can be 
eliminated or means for foreseeing them can be developed. These issues 
will relate to certain criteria of quality, among which is the effect on the 
customer. 

A basic rule is that the customer must not be misled. Some faults, 
whilst not reducing the convenience or reliability of the product, risk 
making a bad impression on the customer and consequently can assume 
great importance, particularly if they are easily detectable. They must 
therefore be avoided, and for this the following scales are adopted. 


Probability of occurrence 


4 Possible jPaales 
3 Improbable 10p. oP 0 
2 Very improbable 10. ee ee 
1 Virtually impossible j aaa (Ni 
Seriousness Detectability 

4 Very critical 4 Very visible 

3 Critical 3 Detectable 

2 Not critical 2 Not very evident 


1 Without effect 1 Undetectable 
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These will be kept in mind by those doing the analysis and are entered 
in the tables of results. Thus the critical issues show up and means for 
preventing their occurrence will appear: all this can contribute to the 
programme of quality improvement. 

We can define a ‘criticality coefficient’ C as 


C = P(probability) x S(seriousness) x D(detectability) 


and the most critical faults will correspond to the highest values of C. 


Example 10 

We consider the application of the FMA method to an adjusting device 
for a headlight beam —- a component that would be supplied to an 
automobile manufacturer by subcontractors. 

There are three items: a nut; an adjusting screw; an adjusting knob. 
Initially the manufacturer sets the adjuster with the vehicle unladen, 
directing the beam correctly with respect to the road (Fig. 2.20). During 
its life, however, the vehicle will be subjected to a variety of loadings 
and the device must therefore perform the function of adjusting the 
beam according to the load. 

Taking C = 24 as the criterion for criticality, the analysis eee D1) 
shows the critical issues to be as in Table 2.8. 


2.2.3 Function analysis 
Function analysis takes account of the relations between 


e@ the need that has to be satisfied and the system being studied, 
together with its functions 
@ the impact of the need on the customer 


Included in a study of reliability it enables 


e the working of the system to be better understood 

@ FMA to be applied 

® communication with the research and development organization to be 
improved 

@ the maintainability and cost controllability parameters to be taken 
into account 

@ failures resulting from links between different components to be 
brought to light in the FMA study 


(a) The method 
The system is considered as a whole and all factors relating to it are 
taken into account. Thus the first things to do are as follows. 
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Figure 2.20 Headlight beam adjustment mechanism. 
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Table 2.8 Critical failure modes and treatment for head- 
light beam adjuster 


Risk Prevention 

Al Breakage of adjusting Strengthen this part 
knob 

<2 Breakage of ring Carry out tests 

3.1 Breakage of angular = Carry out accelerated 
limit stop fatigue and ageing tests 


ee EE EEE EERE ESI 


e Separate the functions from the hardware that realizes them. 
e Identify all the functional components. 


There are three main stages, concerning 


® meeting the need 
@ defining the functions 
© constructing functional block diagrams 


(i) Meeting the need 
The relevant questions are shown in Fig.2.21. This leads to the 
following questions: 


1. Why does this need arise? 
2. How might it be eliminated? 
3. What are the probabilities of the possible ways of achieving (2)? 


In general, (2) will already have been considered by the customer. 


For whom is 
the service 
required? 


With what does 
it provide the 
service? 


Why is this 
action 
needed? 


Figure 2.21 Relationship between various functions of a product or system. 
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(ii) Defining the functions 
A system in a given state is in some kind of contact with its environ- 
ment, to which it provides services of two types. 


1. The services that are its raison d’étre: these are its main functions. 
2. Services that arise in response to the reactions of, or constraints 
applied by, the environment: these are its response functions. 


The main functions correspond to a flow of control across the system 
and are therefore also called flow functions. There are also design 
functions: these are the elementary functions of the components and 
correspond to loops within the system, depending on the design. 

In reliability studies the tasks to be performed are as follows. 


1. Represent the functions and their interrelations by means of a 
functional block diagram. 
2. Quantify the possibilities of breaking the flow, by probabilities. 


(iii) The functional block diagram 
The block diagram is a functional representation of the system in a 
given state of use, showing 


the external environment 

the constituent elements of the system 

the ‘open’ flows, i.e. the main and response flows 

the internal, design, flows 

the contacts, real or virtual, between the system elements 


(b) The different types of function of a product 

Service functions correspond to the needs expressed by the customer; 
they are independent of the technology employed. They must be defined 
with complete clarity in the specification, for on this will depend the 
level of satisfaction of the customer with the product delivered. They 
comprise the following: 


@ the functions that must be provided in order to meet the needs 
expressed; 

@ any additional functions (concerning aesthetics, for example) that are 
provided to improve the customer’s opinion of the product. 


Technical functions are consequential on the service functions and 
depend on the design. In a reliability study an analysis of the technical 
functions enables more to be learned about the consequences of the 
various failure modes. 
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(c) Functional tree, functional block diagram 
The general form of the tree is of the type shown in Fig. 2.22. 


Overall function 
of the product 


1 3 


Figure 2.22 Functional tree for a typical product. 


Figures 2.23 and 2.24 together with Table 2.9 show the application of 
this to Example 10. 
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Controlling the manufacturing 
process 


3.1 VARIABILITY INMANUFACTURED PRODUCTS 


Manufactured products that are supposed to be identical will in fact 
vary, and it is very important to understand this variability; such an 
understanding can lead to 


e reduction in the number of items to be scrapped, 

@ better adaptation of the machines to the production programme, 

® better appreciation of the problems of the production processes in the 
production planning office and 

@ better use of the control charts. 


The study we give here is conducted in terms of an example from 
mechanical engineering (Fig. 3.1); however, the concepts apply equally 
to the manufacture of electronic components and indeed to any type of 
serial production. 

We consider a very simple item, a spacer, manufactured in quantity, 
where a certain dimension is specified; but when samples from a batch 
are measured with a precision gauge the values found are as given in 
Table 3.1. 

These variations can result from a variety of factors: 


temperature changes 

vibration 

positioning in the machine tool 
deformation 

flexure of the machining tool 
wear of the machining tool 


They can be put into two classes: 


1. random variations (RV)—these can be of either sign, and it is not 
possible to predict which. 

2. systematic variations (SV)—the development of these can be pre- 
dicted. 


Fabrication study R&D Office 
Phase: 


Authorised: P. Lyonnet 


Details of operation 
Cf20 = 6+0.2 
routing Cf21 = 16+0.5 


Figure 3.1 Machining instructions. 


Table 3.1 Sample of spacing pieces 


No. Values No. Values No. Values 


5:95 1 5.96 i> “5:93 
25.94 2° 5.96 aw SSGh5 
3 6.00 Bi W503 3 6.06 
4. 5.94 Ato 5008 4 6.01 
5 6.00 5 oe tes Soweisigg 
6 500 6 6.09 6 * 582 
7 307 Id £05.96 6.» 5.82 
Tecwtes 97 TA 5% Md 6.07 
Sy oGt0 eg EP 8 e502 
OF / 5.08 9 5.93 9 ANS OT 
Oe Ae e5.64 10 5.94 10 5.94 
ey  w50Rt Wax, Tencrree x te 5.078 
Woe 016 We O17 W 0.3 


X = mean, W = range. 
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The total variation TV is the sum of these: TV = RV + SV. 

It is important to distinguish between the two types and to measure 
both. If a group of items is taken from a batch, measurements on these 
will show the random variations but will not reveal anything about wear 
or drift. 


3.1.1 Random variations: an example 


After the volume production process has settled down a sample of 20 is 
measured, with results as in Table 3.2. When plotted, these give the 
histogram in Fig. 3.2. The form is characteristic of a random variation; 
the curve of the normal (Gaussian) law is superimposed for comparison. 


Table 3.2 Measured values for 20 samples from a production run 


No. Value No. Value No. Value No. Value 


1 6.10 6 5.80 11 6.20 16 6.00 
2 6.05 7 5.86 2, Pe 17 5.94 
3 5.98 8 5.90 13 6.04 18 6.10 
4 6.00 5 6.01 14 S92 19 5.98 
5 6.00 10 6.05 15 6.01 . 20 ee 


em 


Gauss law 


Number of items 


579 5.85 5.91 5.97 6.03 6.09 Gab 6.21 


Figure 3.2 Histogram showing distribution of values for data in Table 3.2. 
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We shall assume generally that random variations follow the normal 
law: in any particular application this can be checked by applying a 
goodness-of-fit test such as the y? or Kolmogorov-Smirnov test. If the 
parameters of the distribution are the mean m and the standard 
deviation o we can then use the facts that about 95% of the values will 
lie in the range m+2o0 and about 99% in m+3o0. We can find 
estimates for m and o from a sample. The mean for the sample is an 
estimate for m, and if W is the range or spread of values in the sample, 
i.e. the difference between the greatest and least values, an estimate for 
o is W/d,, where n is the number of items in the sample and d, is a 
known function (Table 3.3 gives values of d,,). 


Table 3.3 Estimation of o from range W : 6 = W/d, 


Size of each sample 1/d, ad 
@ 0.886 1.128 
3 0.591 1.693 
4 0.486 2.059 
5) 0.430 2.326 
6 0.395 2.534 
Ti 0.370 2.704 
8 0.351 2.847 
9 0.337 2.970 

10 O25 3.078 
11 0.315 3.173 
ey 0.307 3.258 


Another important property of the normal law is that the mean of a 
sample of size n is itself distributed normally, with mean m and 
standard deviation o/n'”, i.e. the sample mean is the same as the mean 


for a single item but its range of variation is narrower by a factor ee 


3.1.2 Systematic variation: an example 


(a) Introduction 

Ten samples of the same item are measured (a) at the start of a 
production run, (b) after 80 have been made and (c) after 160 have 
been produced. The results are given in Table 3.4. 
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Table 3.4 Samples taken at different stages showing 
systematic change of the mean (x) 


1st sample 2nd sample 3rd sample 
start of after 80 after 160 
series items items 
No Values No. Values No. Values 
1 6.08 1 5.85 1 5.61 
2 5.94 2 Do 2 2 5.68 
3 6.06 8 5.86 3 5.61 
4 6.00 4 5.73 4 Sues) 
5 6.08 5 5.80 5 DOS 
6 5.93 6 5.90 6 5.68 
7 6.09 7 S71 7 Soe 
8 5.92 8 Sefels) 8 5.62 
9 5.98 9 5.2 9 5.51 
10 5.93 10 5.82 10 5.69 
X 6.00 xX 5.79 X 5.61 
Ww Oly W 0.17 WwW 0.18 


Systematic variation is indicated by a steadily changing mean. In this 
example we see that the mean X is decreasing with the number of items 
machined, whilst the range W (the spread about the mean) remains 
effectively constant. This indicates wear of the machine tool, the 
development of which can be represented by a straight line. 

In section 3.2 we shall show that this linear variation of wear can be 
used to construct a control chart for the process. 


(b) An application 
Let TI be the tolerance interval for a measurement of an item; then if 
the variations are random with standard deviation o we must have 

TI = 60 


Otherwise there is a risk that items are rejected. 


Example 1 

For the piece being considered, suppose that the dimension C is given 
by C=5.00+0.20, so that TI=0.40. A sample of 10 gives the 
following measurements: 


5.10, 4.95, 4.93, 5.12, 5:45, 4.90, 5.10, 4.95, 4.85, 5.04 


The range W =5.15 — 4.85 = 0.30; therefore an estimate for o (using 
Table 3.3) is 
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da 3.078. 


This gives 60 = 0.6, which is greater than the tolerance 0.4. Therefore 
we can expect that items will be rejected. 


Example 2 
Suppose that C = 6.00 + 0.025 and TI =0.50. A sample of 10 gives the 
following: 


3.84, 5.A), 6.10, 6.00, 5.97, 5.90, 6.05,°5:98, 6.10, 6205 


w =6.10—5.84=0.26; thus o=0.26/3.078=0.084 and 60=0.5. 
Hence TI = 60, which is a warning that the process should be moni- 
tored. 


Example 3 
C = 4.00 + 0.30 and sample values are as follows: 


4.10, 4.05, 3.95, 4.05, 3.97, 4.03, 3.95, 4.02, 3.96, 4.04 


w = 4.10 — 3.95 = 0.15, o = 0.15/3.078 = 0.048, 60 = 0.29 and TI > 60. 
Therefore there should be no risk of rejects. 


3.2 MONITORING THE MANUFACTURE 


The aim of checking during manufacture is to keep the manufacturing 
process under control and hence to ensure uniform production, whether 
the items manufactured are electronic or mechanical components, food 
products or anything else. To achieve this we must be able to detect any 
deviation from the norm so that we can make the necessary adjustments 
before the process produces items that have to be rejected. 

The parameters monitored can be either attributes — and classed as 
either ‘good’ or ‘bad’ — or properties that can be measured; measure- 
ments of properties give more effective control but are not always 
available. The monitoring uses a graphical presentation called control 
charts; these can be shown to the customer to justify any claim to 
quality of production. 


3.2.1 Control charts using measurements 


The essential requirement is that the parameter used is a measurable 
property, for example length, weight, electrical resistance. As already 
stated, we assume here that the variations observed have a normal, or 
Gaussian, distribution, and so we can use the general results of section 
Bartels 
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When the tolerances on the measurement are known a control chart 
does the following: 


@ it shows up any drift in the measurement concerned, 

@ it enables the intervals at which adjustments should be made to be 
calculated;. 

@ it enables the need for a major resetting to be foreseen; 

© it shows up any increase in the range of variation of the measurement 
concerned, and therefore the need to examine the machine; 

@ jit enables the quality of the manufacture to be assessed. 


If the tolerances are not known the chart provides only the first four of 
these. 

The chart is very much a picture of the manufacture and highlights 
any problems clearly; this is illustrated by the examples in Fig. 3.3. 
When the parameter plotted on the chart is the mean or the median the 
range or the standard deviation should be shown also, to give an overall 
picture of the process. 


(a) Control chart for the mean 
(i) Mean and standard deviation known 

The chart is a plot of the mean values of the parameter we wish to 
control for a series of samples of the same size, n say. We know that if 
the value of this parameter for a single item has a normal distribution 
with mean m and standard deviation o the sample mean will have a 
normal distribution with mean m and standard deviation o/n’?. We 
therefore mark on the chart (Fig. 3.4) 


upper/lower control limits m + 3.090/n 


upper/lower monitoring limits m + 1.960/n¥ 


95% of the points should lie within the monitoring limits and 99.8% 
within the control limits. Thus if a newly plotted point lies on or just 
beyond one of the monitoring limits this is a warning that the process 
needs watching; if it lies on or beyond a control limit some corrective 
action is needed. 


(ii) Mean and standard deviation not known 
When the mean and standard deviation are not known we have to 
estimate these parameters of the law. Suppose that we have r samples 
each of size n and let m;, o; be the mean and standard deviation 
respectively of the values for the ith sample. The estimates are as 
follows: 


xxx xX x 


xx x x xX 


xx x KX x 


X\X xX XX 


wuNK XK & 


xx x x Xx 


xXx xX xX 


ejdwes y}u 


ajdues pig 


g|dwes puz 


a|dwes js 1 


Stable situation; no problem 


‘und uononpold & Ur JoUOD Jo sso] TeoIdk], ¢°¢ aanSry 
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a|dwes puz 
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Drift showing that resetting will soon be necessary 
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Increasing drift and dispersion 
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Upper control limit 


Upper monitoring limit 


x 
A eee ee 
x x x x 
Lower monitoring limit 
Lower control limit 
1 2 3 4 5 6 wh 8 


Sample no. 


Figure 3.4 Control chart: Mean and standard deviation not known. 


=m; 
ite 
; 
1 S03 
= 
Sg 
where b, is given in Table 3.5. If xj, j ='1,2,..., are the values 


measured in sample 7 then 
Pe Fs Ga Seay m;) )1 
es O; ee 
n 


mM; 


n—I1 


The denominator n — 1, instead of n, in the expression for 0; gives what 
is called an unbiased estimate for the sample standard deviation. 

For n < 12 an estimate for 0; is W;/d,, where W, is the range for the 
sample i and d, is the function of n referred to in section 3.1.2(a); 
values are given in Table 3.3. This is quicker, but less accurate. 

When these estimates for m and o have been found they can be used 
as above to construct the control chart. 

Control charts of this kind are used when the tolerance limits are not 
known, or when the spread of values is very narrow. 

If the tolerance interval is known and is large compared with the 
standard deviation, the limits on the chart must be modified to give 
more flexibility in the manufacture. This is discussed in the next section. 


Example 4 
Four samples of a product, each of five items, are measured, with the 
result given in Table 3.6, from which we find 


m; 15.247 15.244 15.238 15.254 


o; 90.018 0.015 0.008 0.011 
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Table 3.5 Estimation of population standard deviation 
from sample value: coefficient b 


Size of each sample _1/b, bn 
2 1.773 0.564 
1.381 0.724 
4 1.253 0.798 
5 1.189 0.841 
6 1.151 0.869 
7 1.126 0.888 
8 1.107 0.903 
9 1.094 0.914 
10 1.083 0.923 
if 1.075 0.930 
12 1.068 0.936 
13 1.063 0.941 
14 1.058 0.945 
15 1.054 0.949 
16 1.050 0.952 1 
17 1.047 0.955 si ene abe 
18 1.044 0.958 | 
19 1.042 0.960 
20 1.040 0.962 e 
/ 
ATG = | (Kee xy] 
21 1.037 0.964 Lae 
1D) 1.035 0.966 
23 1.034 0.967 
24 1.033 0.968 
25 1.031 0.970 
26 1.030 0.971 ee ee 
27 1.029 0.972 Cee ae » On, 
28 1.028 0.973 
29 1.027 0.974 
30 1.026 0.975 


Table 3.6 Measurements of samples from a pro- 
duction run 


Sample 1 Sample 2 Sample 3 Sample 4 


15.250 15.260 15.240 15.240 
157220 15.230 15.230 15.250 
15.240 15.240 15.250 15.270 
15.260 15.260 15.230 15.260 
15.265 {5.250 15.240 15.250 
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The means of these are 15.244 and 0.013 respectively; thus the estimates 
for the population mean and standard deviation are 


_ 0.013 


m = 15.244 o = 0.013 x 1.189 = 0.015 


Hence 
upper/lower control limits = 15.244 + 3.09 x 0.015/V/5 = 15.244 + 0.021 


= 15.265, 15.223 
upper/lower monitor limits = 15.244 + 1.96 x 0.015/V/5 = 15.244 + 0.013 
= 15252 alo.2o3 


from which the control chart can be constructed. 


(b) Modification to take account of known tolerance limits 

If the tolerances are known and the interval is large compared with the 
standard deviation the chart should be modified to take these into 
account, in order to avoid unnecessary adjustments to the process. 

In Fig.3.5, T, and 7, are the upper and lower tolerance limits 
respectively, i.e. the limits between which the measurement must lie; 
lines A are modified control limits and B are modified monitoring 
limits. It can be shown that 


uae 3090 - eA 
n 

T, — B= 3.099 282 2B 
n 


Figure 3.5 Control chart for the mean. 
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If the value of the popuiation standard deviation o is not known, an 
estimate can be found as described above. 


Example 5 
Suppose the tolerance limits are T,, = 122.350 and T, = 122.000 and that 
the estimate for o from samples of five specimens is 0.016. 


3.090 1.960 
= 0.02 : = 
V5 027 3.090 5 


from which it follows that the modified limits are as follows. 
Control: 122.350—0.027 = 122.323 122.000 + 0.027 = 122.027 
Monitor: 122.350—0.035= 122.315 122.000 + 0.035 = 122.035 


3.090 — 


= 0.035 


(c) Control charts for the variations 

The derivation of the control chart for the mean is based on the 
assumption that the variation is stable, and it is therefore necessary to 
ensure that this is so. For this we need to control 


@ the standard deviation (control chart for o) 
@ the spread (control chart for W) 


(i) Control chart for standard deviation 
We suppose that the population standard deviation o is known; when 
this is so it can be shown that, if G is the estimate obtained from a 
sample, the quantity nG/o where n is the size of the sample is 
distributed as y? with n — 1 degrees of freedom. Then if C,, C,, M, and 
M, are the upper and lower limits for control and monitoring respect- 
ively 


Cy een? 
and correpondingly for the other limits, where 
(0.999; n — » ae oe tie 1) 
a le ye 


n n 


Boy 7 


x (0.975; n — 1) 7° (0.025; n — 1) 
Losiale peed esate [Delon “ins det w onal 
n n 
Values of the coefficients B as functions of the sample size n are given 
in Table 3.7. 
If o is not known these limits are given in terms of the estimated 
standard deviation 6 by a set of relations 


C, = BAd etc. 
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and the modified coefficients B’ are given in Table 3.8. 


Table 3.7 Determination of control chart limits: standard 
deviation known 


a acini ni Pe ed We ee 
Size of Ba, Ba 137s Bri 
each sample 


2 0.001 EON 0.022 1.585 
3 0.026 2.146 0.130 1.568 
= 0.078 2.017 0.232 17529 
5 0.135 1922 0.311 1.493 
6 0.187 1.849 0.372 1.462 

7 0.223 1/91 0.420 1.437 
8 0.274 1.744 0.459 1.415 
9 0.309 1.704 0.492 1.396 
10 0.339 1.670 0.520 S79 
11 0.367 1.640 0.543 1.365 
12 0.391 1.614 0.564 1.352 
13 0.413 L59r 0.582 1.340 
14 0.432 1.570 0.598 1.329 
iS 0.450 1552 0.613 1.320 
16 0.467 1535 0.626 1.311 
17 0.482 1.520 0.637 1.303 
18 0.495 1.505 0.648 295, 
19 0.508 1.492 0.658 1.288 
20 0.520 1.480 0.667 1.282 
ZA 0.531 1.469 0.676 1.276 
Za 0.541 1.458 0.684 1.270 
23 0.551 1.449 0.691 1.265 
24 0.560 1.439 0.698 1.260 
25 0.569 1.431 0.704 1.255 
26 0.577 1.423 0.710 1.250 
DH 0.584 1.415 0.716 1.246 
28 0.592 1.408 0.721 1.242 
29 0.599 1.401 0.727 1.238 
30 0.605 1.394 0.731 1.235 


Example 6 
If o(known) = 0.012 and the sample size n = 5, we have 
Cy = :922-<-0:012=-0:028 C, = 0.135 x 0.012 = 0.002 


M, = 1.493 x 0.012 = 0.018 M, = 0.311 x 0.012 = 0.004 
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Table 3.8 Determination of control chart limits: standard 
deviation unknown 


Size of Bu dl Bru ml 
each sample 

2 0.002 4.126 0.039 2.810 

3 0.036 2.964 0.180 2.166 

4 0.098 2.528 0.291 1.916 

5 0.161 2.285 0.370 La 

6 0.215 2.128 0.428 1.682 

7 0.262 ZA0UT 0.473 1.618 

8 0.303 1.931 0.508 1.567 

9 0.338 1.864 0.538 1.527, 
10 0.367 1.809 0.563 1.494 
11 395 1.763 0.584 1.468 
12 0.418 1.724 0.603 1.444 
13 0.439 1.691 0.618 1.424 
14 0.457 1.661 0.633 1.406 
15 0.474 1.635 0.646 1.391 
16 0.491 1.612 0.658 13377 
7 0.505 1.592 0.667 1.364 
18 0.517 {571 0.676 19852 
19 0.529 1.554 0.685 1.342 
20 0.541 1.538 0.693 1.333 
21 0.551 1.524 0.701 1.324 
Ds 0.560 1.509 0.708 1315 
23 0.570 1.498 0.715 1.308 
24 0.579 1.487 0.721 1.302 
25 0.587 1.475 0.726 1.294 
26 0.594 1.465 0.731 1.287 
Pap 0.601 1.456 0.737 1.282 
28 0.608 1/447 0.741 ft 76 
29 0.615 1.438 0.746 1271 
30 0.621 1.430 0.750 1.267 


An important point to note here is that if the sample standard deviation 
falls to one of the lower limits, or below, this is not a signal to intervene 
in the manufacturing process but rather to inspect the measuring 
equipment so as to ensure that it is not giving faulty readings. A real fall 
in the standard deviation — i.e. in the spread — is relatively rare. 


(ii) Control chart for range 
There are corresponding coefficients D, D’' for the cases when the 
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population standard deviation o is and is not known respectively. When 
o is known 

C= Bro etc: 
and when a is not known 

C, =Daw ete. 


where W is the sample range or the mean of the ranges of a number of 
samples. Values of D, D’ are given in Tables 3.9 and 3.10. 


Table 3.9 Control chart for range: standard deviation 


known 
Size of De Dg Dx Dm 

each sample 
2 0.00 4.65 0.04 Sol7 
3 0.06 5.06 0.30 3.68 
4 0.20 Sol 0.59 3.98 
5 0.37 5.48 0.85 4.20 
6 0.54 5.62 1.06 4.36 
i) 0.69 S773 175 4.49 
8 0.83 5.82 1.41 4.61 
9 0.96 5.90 1355 4.70 
10 1.08 5.97 1.67 . 4,79 
11 1.20 6.04 1.78 4.86 
12 1.30 6.09 1.88 4.92 

Table 3.10 Control chart for range: standard deviation un- 

known 

Size of Dos a Dine td 

each sample 
Z 0.00 4.12 0.04 2.81 
3 0.04 29 0.18 Zo, 
4 0.10 2.58 0.29 193 
©) 0.16 2.36 0.37 1.81 
6 0.21 Q22 0.42 D372, 
ie 0.26 pam) Pe 0.46 1.66 
8 0.29 2.04 0.50 1.62 
9 O32 1.99 0.52 1.58 
10 0.35 1.94 0.54 1.56 
Ld. 0.38 1.90 0.56 1:53 
iz 0.40 1.87 0.58 Lok 
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Example 7 

If the range W = 0.122 for a sample of size n = 10 the four coefficients 
D., etc. are 1.94, 0.35, 1.56, 0.54; multiplying the observed range 0.122 
by these gives the limits 


control (0.234, 0.043) 
monitor (0.195, 0.065) 


For practical use of the chart it is convenient to shade the bands 
between the control and monitor limits, as in the example in Fig. 3.6. In 
the case of the chart for the standard deviation only the upper limits 
need be shown. 


Upper control limit 


ON 


Lower monitoring limit 


Sample no. 
a 


(a) 


iy control limit 


Figure 3.6 (a) Control chart for the mean (example 7). (b) Control chart for the 
standard deviation (example 7). 
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(d) Control chart with directly plotted individual values 

With this form no calculation at all is done, the individual values being 
plotted directly (Fig. 3.7). The rule here is that the occurrence either of 
two successive values in a shaded band or of one value outside either 
control limit is a signal to take action. 


Upper control limit 


Upper monitoring limit 


x x 
x x x 
x x x x = is x 
2 x x x x x . x 
Lower monitoring limit x 
Lower control limit 
1 2 / 3 4 5 6 fi 8 


Sample no. 
ee 


Figure 3.7 Control chart for individual values. 


In constructing this form of chart the assumption is made that the 
tolerance interval TI is six times the population standard deviation: 
TI = 60. Thus the estimate for o is 6 = TI/6. Further, the distances d, 
and d, of Fig.3.8 can be found using d; = P; x TI and d, =P, x TI 


where P, and P, are found from Table 3.11. 


Upper control limit d, 
Upper monitoring limit 


Lower monitoring limit 


Lower control limit 


T 


Figure 3.8 Control chart for the number of defects allowed. 
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Table 3.11 Determination of limits for Figure 3.8 


ey P = 0.5% P = 0.27% 
n (U=2.6) (U = 2.8) (U = 3) 
3 P, 2.6% 6.5% 9.3% 
P, —-:16.3% 19.1% 21.1% 
4 P, + 0.7% 4.7% 7.6% 
P, 13.4% 16.4% 18.5% 
5 P,  — 0.8% 3.4% 6.4% 
Boy gale 14.4% 16.7% 
6 Pi = 1.9% 2.3% 5.4% 
Py cee 6% 12.9% 15.3% 
7 P,  — 3.0% 1.4% 4.5% 
Pr 8.3% 11.7% 14.2% 
8 P,; — 3.8% 0.6% 3.8% 
Py © 72% 10.7% 13.2% 
9 Py = 4.5%  — 0.1% 3.2% 
P, «63% 9.9% 12.4% 
10 Py = 5.1%  — — 0.6% 2.6% 
P; gcd 4% 9.1% 11.7% 


Example 8 
Find the limits for the individual value chart with 


T=) 1 50),250 T, = 150.000 Te = Os200 
n=.6 


and the percentage P of defectives acceptable equal to 0.5. 


For these values Table 3.11 gives d; = 0.023TI, d, = 0.129TI; hence 


C,, = 150.250 — 0.023 x 0.250 = 150.244 
M,, = 150.250 — 0.129 x 0.250 = 150.218 
M, = 150.000 + 0.129 x 0.250 = 150.032 
C, = 150.000 + 0.023 x 0.250 = 150.006 


3.2.2 Control charts for attributes 


Not all the properties that affect quality can be measured: this is 
especially the case where any kind of aesthetic judgement is involved, as 
for example with the coachwork of an automobile. The control criteria 
must then be qualitative and expressed in terms such as ‘good/bad’, 
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although of course some quantitative control may be involved, such as 
some dimension being or not being within the allowable limits. In these 
circumstances it is useful to distinguish between control charts for 


e@ the number of defective items 
@ the proportion of defective items 
@ the number of defects per item 


(a) Control chart for the number of defectives in a sample, the 
proportion in the population being known 

The control and monitor limits are defined as those within which 99.8% 
and 95% respectively of the population lie; their values are calculated 
on the assumption that the fraction P of defectives in the population is 
known. 

If the sample is of n items and the total population is N items then 
(see Chapter 6) if n/N <0.1 the binomial law can be used to calculate 
the probability of any number of defectives in the sample; otherwise the 
hypergeometric law must be used. 

Assuming that the binomial law is applicable, the upper control limit 
C, is defined as follows (we are not, of course, interested in a lower 
limit): 

Cy 
S ) Pk(1 — P)"-* = 0.998 
k=0 ; 


The upper monitor limit is given by 


M, 
> ) PK(1 — P)"-* = 0.975 
k=0 


These limits are now integers as we are dealing with discrete variables; 
consequently we may not be able to achieve the probability value 0.998 
and 0.975 exactly. 

Values for the limits can be found by using Tables 3.12 and 3.13 
which have been computed from the binomial and Poisson distributions. 


Example 9 
The proportion of defectives is stable at 3%; we want to know the 
control and monitor limits for a sample size of 15. 


Looking down the columns headed n = 15 of Table 3.13 we see that the 
nearest we can get to 3% is with 


C, = 3, corresponding to P = 3.1% 
M, = 1, corresponding to P = 1.9%(M,, = 2 corresponds to 4.3%) 
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(b) Control chart for the proportion of defectives 

A control chart for the proportion of defectives is preferable when the 
size of the sample can vary. For a sample of n the corresponding limits 
are 


C= 


n n 


and if as before k is the number of defectives in the sample the quantity 
to be controlled is 


Dp =k jn 
Thus for the above example, Ci, = 3/15 = 0.2, M/, = 1/15 = 0.067. 


3.3 INTERVAL BETWEEN CONTROL ACTIONS 


One of the problems facing the quality engineer is to predict the 
intervals at which control actions should be taken; this will depend on 
the way in which the equipment can get out of adjustment, which can be 


@ by random changes or 
@ by systematic drifts, due to wear, overheating etc. 


3.3.1 Random variations 


A study of the production statistics will enable the desirable interval 
between control actions to be estimated; the control charts will indicate 
when the process is going out of control. However, it is important to 
find the reasons for these changes, and a study of the correlations 
between the changes and the parameters of the process will help here. 


3.3.2 Drifts 


This type of change is easier to correct because, given enough observa- 
tions, regression analysis (linear, polynomial, exponential or loga- 
rithmic) can be used to model the law describing the change. 

Tool wear is a cause of drift in mechanical engineering production, 
and tribological studies have shown that this is practically linear over the 
active life of the tool: this is illustrated in Fig. 3.9. Thus regression can 
be used to model the change. 

Tool life can be expressed as a distance, the distance travelled either 
by the tool itself or the metal or other material that is being machined. 
If y is the amount of wear (also measured as a distance) when the tool 
life is x, the relation assumed _is 


y=ax+b 
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Wear in 1/100(mm) 


Burn-in Active life Old age Distance 
travelled by the 
tool in metres 


Figure 3.9 Typical life history for machine tool wear. 


If we have n pairs of observations (x;, y;), with means X and Y 
respectively, the least squares fitting of the straight line gives (Fig. 3.10) 


= Xi; —= nXY 
ESS 

= (x)? — nx? 
b=Y-aX 


This can now be used to predict the stage at which the wear will have 
reached the maximum that can be allowed. 


Upper control limit 


Sample no. 


Figure 3.10 Control chart with successive samplings showing cumulative effect of 
wear. 
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Example 10 
Find the regression line for the following record of values, and the tool 
life (in metres) if the allowable wear is 2 mm. 


Here y (in millimetres) is the relevant dimension of the piece machined 
when the tool life is x (in metres). 


yy 35:10' 35.15 3516 35,17535.19"35.20°35.29 Si 24035 30835240 
x. "00 aio. BUG 0Sia 2900". SET “7 Sid*5 2515 ~ 520° 53055 


Using the results given above we find the line to be 
y = 30.162 + 9.877(x/1000) 
Thus a change Ax results in a change Ay where 
Ay = 9.877(Ax/1000) 
and so, for Ay = 2, Ax = 2000/9.877 = 202 m. 
Example 11 involving the control chart 
We wish to establish a control routine, with chart, for an item to which 


the following apply: tolerances T, = 75.50, T, = 76.50; sample size 
n = 5, Measurements on a pre-production sample give 


1290 75.93,. 19.90 73.92 “TSOBR AS. 95% FS 15.94 
TOD T3290 T9502 NS. 90" 1S. IF IG OFS. I9 


Measurements made in order to determine tool wear are given in Table 
3.14, where y; is the dimension for the ith piece machined. The time 
required to machine one piece is 12 min. 


Table 3.14 Data for example 11 


L Yi u Ji 

1 75.000 42 75.042 

3 75.003 80 75.080 
10 75.009 90 75.081 


ZS 75.026 100 75.101 
30 75.031 150 75.149 


We first estimate the population standard deviation o: using the formula 
given in section 3.2.1(a) with the values for the pre-production sample 
we find 
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6 = 0.0294 


Next, we find the limits for the control chart for the mean by using the 
method of section 3.2.1(b): 


= = 0.0489 = T|A Die 2 a ae — = 0.062 = T, 
The chart limits are therefore 

C,= T, — T,A = 76.500 — 0.049 = 76.451 

M, = T, — T,B = 76.500 — 0.062 = 76.438 

M,= T, + T,B = 75.500 + 0.062 = 75.562 

C,= T, + T,A = 75.500 + 0.049 = 75.549 


T,A =36- 


Next we find the chart limits for o, for which only the upper values are 
of interest: 


C, = Bud = 1.922 x 0.0294 = 0.056 
M,, = Bm = 1.493 x 0.0294 = 0.044 


(from Table 3.7). We can now construct the control charts. 
To determine the frequency of control action, we first find the 
amount of wear that can be allowed: 


M = TI — 60 = (76.500 — 75.500) — 6 x 0.0294 = 0.824 
The regression line for the wear data is found to be (Fig. 3.11) 
y = 9.85 x 107* + 74.999 


and so the number of pieces machined when the wear has reached 0.824 
is 0.824/9.85 x 107* = 836. 

Finally, since the time required is 12 min per piece, the tool should be 
reset at intervals of 836 x 12/60 h = 167 h of working life. 


10 REM REGRESSION LINE 


20 CLS 

30 LOCATE 8,15: INPUT ‘How many (X,Y) pairs do you wish to consider”; 
N 

40 CLS a 

50 LOCATE 4;51: PRINT “X pe 

60 DIM X(N): DIM Y(N) 

70 CLS 

80 FOR I=1 TON 

90. LOCATE I°2, 5: PRINT “PAIR” :1:: INPUT “Value. ef X=’, X(1) 


Figure 3.11 BASIC program for finding regression lines: only X;, Y; pairs need 
be entered. 
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100 LOCATE I*2,40: INPUT “Value of Y= _ __”, Y(I) 

110 NEXT I 

120) CLES 

130 LOCATE 4,6: PRINT “Number of pairs (X,Y):”; PRINT N 
140 LOCATE 4,51: PRINT “X Nae 


150 FOR I=1 TON 

160 LOCATE 5+1,50: PRINT X(1), Y(D) 

170 LOCATE 5+1,45: PRINT CHR$(179) 

180 PRINT: X=X+X(D: Y=Y+Y(I): X2=XK2+X(I) 2: 
Y2=Y2+Y(I) 2: XY=XY+X(I)*Y(1) 

190 NEXT I 

200 LOCATE 6,1 

210 PRINT “Sum of Xs D8 

220 PRINT “Sum of Ys ag 

230 MX=X/N: MY=Y/N 

240 PRINT‘ Mean (X) os VEX. 

250 PRINT “ Mean (Y) 725 MEY 

260 Print ‘Sum of squares (X) cig. © 

270 Print “Sum of squares (Y) Bae 

280 MXY=XY/N 

290 VX=X2/N—-MX 2: VY=Y2/N-MY 2 


300 PRINT “Variance X TEEN DG 

310 PRINT “‘Variance Y tae VY. 

320 COV=XY/N-(MX*MY) 

330 PRINT ‘Covariance (X,Y) SOON 


340 RO=COV/(SOR(VX)*SOR(VY)) 

350 PRINT USING “ Correlation coefficient #.###’;RO 
360 PRINT “Regression line”’ 

370 PRINT “ ---------- i 

380 A=COV/VX: B=MY—R*MX 

390 PRINT USING “ Coefficient A : #.4##”; A 
400 PRINT USING “ Coefficient B : #.4##”; B 
410 FOR I=1 TO 19 

420 LOCATE 2+I,1: PRINT CHR$(186) 

430 LOCATE 2+1,71: PRINT CHR$(186) 

440 NEXT I 

450 LOCATE 22,1: PRINT CHR$(200) 

460 LOCATE 2,1: PRINT CHR§(201) 

470 FOR I=1 TO 70 

480 LOCATE 22,I+1: PRINT CHR$(205) 

490 LOCATE 2,I+1: PRINT CHR$(205) 

500 NEXT I 

510 LOCATE 2,71: PRINT CHR$(187) 

520 LOCATE 22,71: PRINT CHR$(188) 

530 LOCATE 23,10: INPUT “Any more? (Y/N); A$ 
540 IF A$=Y THEN 20 ELSE 550 

sI0 -END 

Figure 3.11 (continued) 
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Quality control of goods 
recelved 


Quality control of goods received is an important component of any 
quality programme; its purpose is to filter out any below-standard items 
or materials delivered to a part of the enterprise 


e from outside, by any of its suppliers, or 
e@ from inside, by another unit or department. 


The ways in which control is usually exercised are 


@ by attribute(s): the item concerned is rated ‘good’ or ‘bad’ 
according to some criterion, and as with 
control charts the decision on whether to 
accept or reject the delivery is based on 
the number of ‘bad’ items in the sample. 

® by number of defects per the mean number is found, by sampling, 


item: and the decision to accept or reject is 
based on this. 
e@ by measurements: the property on which the decision is to be 


based must be measurable, as when used 
for a control chart; the sample mean and 
standard deviation can be used as criteria. 


4.1 CONTROL BY ATTRIBUTES 


4.1.1 Types of sampling 


(a) Simple sampling 
A single sample is taken and the decision whether to accept or reject the 
batch is based on the number of defectives found. 
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Sample of n items 


‘ 
| ie sa ee 


No. of defectives < A No. of defectives = R (= A+1) 


Accept batch Reject batch 


Here A is the acceptance criterion, R the rejection criterion. 

This is the simplest control procedure to implement; however, it is not 
optimal from an economic point of view, and for batches that are 
definitely good or definitely bad the decision can be based on a smaller 
sample, thus reducing the cost. This can be achieved by double or 
multiple sampling. 


(b) Double sampling 

A first sample is taken as before; depending on the result either a 
decision is reached or a second sample is taken and the decision is based 
on the result of the two combined. 


Samples of n, items, k, defectives found 


| 


Accept | Reject 


Sample of n, items, k, defectives in (n,; + n) 


k, = A, k, 2 R, (= A, + 1) 
Accept Reject 


(c) Multiple sampling 
This process can be repeated: the criteria for the second sample are 


k, < A, accept 
k, = R, reject (R, # A, + 1) 
A,< k,<R, take a third sample n3: 
k; defectives in n, + n,+ 73 
kz; = Az; . ~ accept 
k, = R; reject (R; = A; + 1) 
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and so on, up to eight samples (this number is defined by AFNOR). 
The values of the criteria A,, R, etc. are found by using the standard 
statistical distributions — binomial, hypergeometric, Poisson. 


4.1.2 Laying down a control procedure 


The calculations to define any of a number of control procedures have 
been done and the results are given in the French standard NF 022 X; 
this greatly simplifies the task of implementing the procedure in any 
particular case. Alternatively, a software package can be used. 

The following have to be specified for the procedure: 


e the type of control to be used (by attributes or by measured 
properties) 

@ the sample size 

e@ the method of sampling 

@ the relation between the result of the sampling and the decision taken 


The type of control will depend on whether the property of interest is 
or is not measurable. ; 

The sample size will depend on the size of the batch and on the 
stringency with which the control is to be applied; Tables 4.1 and 4.2 
give standard recommendations, in which 


@ levels $1, S2, S3 and S4 apply to military applications only, 
@ level I is for relatively relaxed control, Level II is normal and Level 
III is for relatively strict control. 


Through the code letter given in Table 4.1, Table 4.2 gives the values of 
n, for simple sampling, n, for double and n; for higher multiplicities. 
The method of sampling will have to be agreed by the two parties. 
For each sampling scheme the decision criteria are functions of the 
acceptable quality level (AQL), expressed as the fraction or percentage 
of defective items in the batch (see Fig. 4.5 later). 


4.1.3 Risks borne by the supplier and by the customer 


It is important to realize that there are inherent risks of making the 
wrong decision in any control system based on sampling, and these must 
be accepted by the two parties. 

Suppose the true (unknown) fraction of defectives in the batch is p, 
and that the batch is acceptable if p < p, and not acceptable if p = pp. 
If k is the number of defectives found in a sample of n items, then with 
the notation of section 4.1.1 this means A = np,, R = np». The risks are 
as follows: 
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100% 


0 p% limit of defects 


Figure 4.3 Idealized control system with 100% efficiency. 
100 


80 


60 


Ov 


40 >» 


20 


100 


50 100 150 200 250 300 350 400 450 500 
(b) 


Figure 4.4 Efficiency curves for simple sampling plans (letter code E in Table 
4.2a). 


The probability p, of accepting a batch containing a fraction p of 
defectives is the probability that the number k of defectives found in a 
random sample of n items does not exceed a stated number A; this is 
expressed formally as 
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4.1.4 Changing from one level of control to another 


Assuming that the process starts in normal level, the rules are: 


if normal and 10 successive batches are accepted change to relaxed 
2 successive batches are rejected change to strict 

if relaxed and 1 batch is rejected change to normal 

if strict and 5 successive batches are accepted change to normal 


This is shown diagramatically in Fig. 4.2. 


|: relaxed control 2) 


||: normal control 
Ill: strict control 


Figure 4.2 State change diagram for control levels. 


4.1.5 Efficiency curves 


A control procedure can be regarded as a filter that allows only those 
batches having a percentage of defectives not exceeding a certain value 
P, to pass through. A perfect, or 100% efficient, procedure would have 
the characteristics of Fig. 4.3, passing all batches having percentage 
defective P < Py and rejecting all with P > Py. As we have seen, this is 
not attainable if the procedure is based on sampling; the curve of the 
fraction p, accepted (or percentage P, = 100p,), based on one or other 
of the standard statistical laws, has the general form of Fig. 4.4. The 
closer this curve approximates to the rectangular form of Fig. 4.3 the 
more efficient is the procedure; but in general this increase in efficiency 
is gained at the cost of increasing the sample size n. 
The statistical laws used in determining the curve are as follows: 


© hypergeometric when the batch N is small and sampling is without 
replacement; 

@ binomial when n/N <0.1, where n is the sample size; 

® Poisson when the batch size is large but unknown. 
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100% 


p% limit of defects 


Figure 4.3 Idealized control system with 100% efficiency. 
100 

80 

60 

40 


20 


50m O08 15029200) 6250; SS005350! 400 450 500 
(b) 


Figure 4.4 Efficiency curves for simple sampling plans (letter code E in Table 
4.2a). 


The probability p, of accepting a batch containing a fraction p of 
defectives is the probability that the number k of defectives found in a 
random sample of n items does not exceed a stated number 4A; this is 
expressed formally as 
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Pa = prob(k < Aln, p) 
Let p, be the greatest fraction of defectives that is acceptable; then p, 
is given in terms of A and p, by the following: 


1. For n/N > 0.1: hypergeometric law 


_"S (Npx) (NG-pd) |(N 
Pav ee\ ke n—k n 
2. For n/N <0.1: binomial law 


k=A 
n eae 
Da = > (7) pia mm P1) ‘ 
k=0 
3. For N large, p < 0.1: Poisson law 

k=A 


exp(—m)m* 
Pa = SS fal 
k=0 : 


where if is the sample size, m = np. 

The efficiency curve can also be found experimentally, by testing 
numbers of batches with different but known fractions of defectives. 
Then 
no. of batches accepted 


Parle, Sere niches ceeted 


for each value of p. Efficiency curves relevant to current practice are 
given in the Standards documents; Fig. 4.4 and Table 4.3 are examples. 


Table 4.3 Levels of acceptable quality 


Level of acceptable quality (normal control) 


Ps 1.0 4.0 6.5 10 
p (percent defective) 
99.0 0.077 1.19 3.63 70 
95.0 0.394 2.81 6.83 11.3 
90.0 0.807 4.16 8.80 14.2 
75.0 2.19 7.41 13.4 19.9 
50.0 5:19 12.6 20.0 21:9 
25.0 10.1 19.4 28.0 36.2 
10.0 16.2 26.8 36.0 44.4 
20 20.6 31.6 41.0 49.5 
1.0 2928 41.5 50.6 58.7 
is) 6.5 10.0 


Level of acceptable quality (reinforced control) 
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Example 1 

Batch size = 80, AQL = 1%. 

From Table 4.1 the code letter for this batch size and a normal level of 
control is E (column IJ); then from Table 4.2(a) the sample size is 13 
and the criteria are 


(accept) A = 0 (reject) R =1 


This specifies the control procedure; the efficiency curve for 1% AQL is 
given in Fig. 4.4(a). 


4.1.6 Quality improvement resulting from control 


The control procedure can be operated in such a way that the stream of 
outgoing items (i.e. after the procedure has been applied) has a higher 
quality than the incoming stream. That is, if p and p’ are the fractions 
of defectives before and after control, then p’ < p. This can be achieved 
by testing all the rejected items and adding those found to be free of 
defects to the output stream. 

Let p, be the fraction of batches that are accepted by the control 
procedure: as we saw in section 4.1.5 this is a function of the (unknown) 
fraction p of defective items in the batch. Thus if we test r batches of N 
items each we accept p,r batches and reject (1— p,)r. The output 
stream at this stage consists of p,rN items among which are pp,rN 
defectives. 

The reject stream consists of (1 — p,)rN items, of which p(1 — p,)rN 
are defective and (1 — p)(1— p,)rN are free of defects. If we remove 
all the defectives and add the defect-free items to the ‘accept’ stream we 
have a total outgoing stream of p,rN+(1-— p)(1—p,)rN = 
(1— p+ pp,)rN items, among which are pp, defectives. 

Thus the outgoing quality, the fraction p’ of defectives in the output 
stream, is 


at. PP a 

1— p+ ppa 
since 1 — p > 0. If, as expected, p is small and p, is close to unity, p’ is 
close to pp. 

Since p, is a function of p and can be calculated for a range of values 
of p when the details of the procedure are known (see section 4.1.5), p’ 
also can be calculated as a function of p. The expression shows that 
p' = 0 when p = 0, which is obviously true (all the items are good), and 
p' =p, when p =1; but then p, = 0 (all items are bad and no batches 
are accepted), so p’ = 0 again: not a meaningful result, since there are 
now no items in the output stream. It follows, however, that p’ will 


Pp edd 
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have a maximum value for some value of p between 0 and 1 (i.e. of the 
percentage P between 0 and 100), which corresponds to a minimum 
outgoing quality after the test. 

Figure 4.5 gives the curve of p’ as a function of p for a control 
procedure of simple sampling with sample size n = 20 and acceptance 
level A = 3 defectives. 


Average quality after control procedure 


4 5 6 7 8 9 
Percentage of defectives in batch 


Figure 4.5 Effect of control: outgoing quality. 


4.1.7 Average number of items tested 


The average total number of items tested can be computed as a function 
of the true quality of the batches, for simple, double and multiple 
sampling procedures. This enables an economic choice of procedure to 
be made. 

For simple sampling the size is n, the value chosen. 

For double sampling let p; be the probability of needing to take a 
second sample. 


k=R, 
pi = prob(A,; < k, <R)) = > 2) Dl pe 
kA, 


where p is the true fraction of defectives in the batch. Then if n, and n, 
are the sizes of the first and second samples respectively the average 
number tested in this procedure is 

Ny = Ny + P aN 


For multiple sampling with an obvious extension of the notation the 
result is 


Sequential testing: Wald’s test 95 
Nay = Ny + pany + pny + os. peng 


If the procedure includes sorting the rejects to remove the defective 
items the average numbers handled become 


e for simple sampling 


Haz Sap all pd Fi p,)n 


where p, is the probability of acceptance; 
e@ for double sampling 


one = Pan i Palny a N) 17 (1 = [Dat = Pa)N 


where p.1, Pa. are the probabilities of acceptance for the samples ny, 
ny respectively, and N is the number of items checked. 


4.2 SEQUENTIAL TESTING: WALD’S TEST 


Wald’s test was developed with the aim of reducing the number of items 
tested; the principle is that items are drawn and tested one after another 
and the decision to accept or reject the batch is taken when one or other 
of two conditions is satisfied. 

The procedure is shown graphically in Fig. 4.6. Items are tested in the 
order 1,2,3,... and the accumulated number of defectives after each 
test is plotted. This gives a stepped line, moving always either to the 
right or upwards, and the decision is taken when the line meets one or 
other of the ‘accept’ or ‘reject’ boundaries. 


No. of defectives 


1 2 See 5 eo) 7 Se 9) 10" Hl 
No. of items tested 


Figure 4.6 Limits for sequential sampling (Wald test). 
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4.2.1 Determining the boundaries 


The boundaries are determined by the two risks w and £ to the supplier 
and customer respectively (cf. section 4.1.3) and the corresponding 
probabilities p,, p2. The equations of the lines are as follows: 


(accept) L1 a, = h, + Sn 
(reject) L2 lr, = hy + Sn 


where 
S= K tos(; -- ad 
hi= Ktos(*] 
h,= K tog(* =f) 
Bio [Paley ea eel 


and logarithms are to base 10. 


4.2.2 Efficiency 


An estimate of efficiency is given by the curve shown earlier as Fig. 4.1. 


4.2.3 Scoring procedure for Wald’s test 


Instead of the graphical method just described a scoring procedure can 
be used, as follows. As before, items are tested one after another. 
We define a score H and two critical values H,, H): 
math h, 


LS S eS 


with h,, h, and S as before. Initially H = H, (the ‘handicap’); after 
each successive test 


add 1 to the current value if the item is ‘good’ 
subtract (1 — $)/S if the item is ‘bad’ 


and accept when H = H, or reject when H = 0. 
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4.3 CONTROL BY MEASURED PROPERTIES 


For control by measured properties, obviously the property chosen as a 
criterion must be measurable; in general this method needs a smaller 
sample size than testing by attributes, for the same efficiency. 

We make the basic assumption that the values of the measurement 
chosen are distributed normally (Gaussian) with mean m and standard 
deviation 0; we assume also that the true values of m and_o for the 
batch are not known and have therefore to be estimated from the 


sample. 
With a sample of values (x,, x2, ..., X,) the estimates for m and o 
are 
aR Bi 
i= — 
n 


ae Ay? i. 
B-Newiar ee ae 
eal! 
Let T,, be the upper tolerance, i.e. the maximum acceptable value for 
the measure. Then the fraction p of defective items in the batch is given 


by the probability that a value exceeds T)\, i.e. 


= rae tll 
P= pro 6 


where, as usual, U is the reduced normal variate (x — m)/o (Fig. 4.7). 
p is related directly to the AQL and is easily found from the table for 
the Gaussian distribution. 

The calculations involved here are the same as for the control charts 
based on measured properties (cf. Chapter 3, section 3.2.1). 


| 
I 
m A 


Figure 4.7 Upper tolerance level (Gaussian distribution). 


4.4 SAMPLING PROCEDURES 


There are many different ways in which samples can be taken, and the 
appropriate one to use in any particular case will depend on how the 
batch has been assembled: the sample should always reflect the reality 
of the situation as closely as possible. 
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In random sampling items are chosen with the help of random 
numbers, which can be obtained from tables (Appendix 8 gives a short 
table) or generated by a computer program — standard programs for this 
are available. This method is used when it can be assumed that the 
batch is homogeneous. 

In stratified sampling the items in the batch come from several 
different sources (‘strata’) and it is desirable to take this into account in 
the sampling. This can be done by drawing random samples from each 
source and combining them to form the sample for the test. Examples of 
this situation are batches consisting of items made in different factories 
or in different runs in the same factory. 

Two-level sampling consists in first selecting large samples from the 
primary units in the batch and then constituting the sample at random. 


5 
Cause-and-effect analysis 


5.1 THE ISHIKAWA CAUSE-EFFECT DIAGRAM 


5.1.1 General principles 


When a manufacturing process is being monitored the first sign that 
something is wrong is the production of items that have to be scrapped; 
if further unwelcome effects are to be avoided the real causes of this 
must be discovered. 

Since manufacturing processes often use complex systems, as many 
people as possible who are able to contribute to solving the problem 
should participate in the investigation, and in particular the users of the 
system. When the group has got together, notice should be taken of all 
suggestions concerning the loss of quality: 


variability of the raw materials; 

variability of the machines involved; 

changes in the workforce; 

changes in the working environment — e.g. from day shift to night 
shift; 

changes in working practices. 


5.1.2 Application of these to the problem 


Experience has shown that the causes of any effect in a manufacturing 
enterprise can be grouped into five main classes, which can be repre- 
sented as a basic cause-and-effect diagram (Fig. 5.1). In the investigation 
each possible cause is recorded on this diagram. The next step is to 
establish the validity of the assertions and the relative importance of the 
various possible causes. They cannot all be investigated at once and so 
they must be put in order, and for this a scheme of weighted voting is 
helpful: each participant gives a weight to each cause and the causes 
that receive the greatest total weights are studied first. 
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Materials 


Machinery 


Figure 5.1 Typical cause/effect diagram for manufacture. 


Example 1 
The case studied is chosen to show how very general this approach is to 
the problems of industry. The problem is 


‘the coffee machine is giving bad coffee’ 
The possible causes listed are as follows: 


water is cold; 

water is too hot; 

water is polluted; 

water is chlorinated; 

poor brand of coffee; 

wrong amount of coffee; 

coffee badly ground; 

poor quality sugar (cane, white); 

too much sugar; 

oil on beaker; 

plastic beaker; 

water temperature; 

wrong amounts/proportions of ingredients; 
machine needs attention; 

machine inconvenient/unattractive; 
location of machine wrong — noise, dust etc. 


Entering these on the basic diagram of Fig. 5.1 we get Fig. 5.2. 
There were five participants in the study; each weighted each possible 
cause on a scale of 0-20 and the highest total scores were as follows: 


@ water quality (hard or soft) (80) 
@ brand and quantity of coffee (65) 
@ beaker (50) 
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Convenience Coffee type 


Appearance 


Water 


Portion control 


Coffee brand 


State of adjustment Sugar 


Water temperature 


Bad coffee 


Location 


Maintenance 


Coffee (ground, instant, liquid) 


Methods 


Figure 5.2 Cause/effect diagram for the coffee machine problem. 


Having agreed on a short list of possible causes the usual next course is 
to make tests and experiments and carry out statistical analyses such as 
correlations, analysis of variance etc. 


5.2 PARETO OR ABC ANALYSIS 


Pareto or ABC analysis is an investigatory tool that enables the quality 
assurance service to assign priorities to the possible sources of quality 
defects — examination of rejects, for example, is the most expensive 
action. It can also be used to assess the improvement that has been 
achieved in any process by comparing the ABC curves for different 
dates. We describe it here in the context of quality in manufacturing 
production, and thus of rejects, but it is of much more general 
application. 


5.2.1 The method 


Rejects are put into classes according to some criterion and the classes 
are arranged in decreasing order of the costs they incur. This is then 
represented graphically — the ABC graph — with accumulated percen- 
tage of costs plotted against the accumulated percentage of types. Figure 
5.3 shows the typical form of this curve. 
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Figure 5.3 Pareto (ABC) analysis. 


Zone A is the top priority zone. In most cases it is found that about 
20% of the rejects account for about 80% of the costs. 

Zone B contains the next 30% of rejects, which account for about 
15% of the costs. 

Zone C: the remaining 50% of rejects account for the remaining 5% 
of the costs. 


Example 2 

The costs associated with the rejects produced by the various machines 
in a manufacturing unit were as shown in Table 5.1. Putting these in 
order of cost we obtain Table 5.2 which, plotted, gives the curve of 
Fig. 5.4. The conclusion is that machines 11, 10, 1, 8, 9 and 3 should be 
examined first; restoring these to proper adjustment would save nearly 
80% of the costs due to rejects. 


Table 5.1 Costs per machine 


Machineno. Costs Machineno. Costs 


(£100) (£100) 
1 100 8 80 
2 32 9 55 
3 50 10 150 
4 19 11 160 
5 4 12 5 
6 30 13 10 
7 40 14 20 


V——————————————————————— 
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Table 5.2 Data for the production of the curve 
in figure 5.4 


Machine Costs Accumulated Percentage 


11 160 160 2A.2 
10 150 310 41.0 
1 100 410 54.3 
8 80 490 64.9 
S) 55 545 ges 
3 50 595 78.8 
7 40 635 84.0 
2 32 667 88.0 
6 30 697 92.0 
14 20 rae 95.0 
4 19 736 hes 
13 10 746 98.8 
AD 5 ip OS 
3) 4 F355 100.0 


0 23.6 50 100% 


Figure 5.4 ABC analysis for Example 2. 


Experience has shown that in general the ABC curve will have one or 
other of the three forms of Fig. 5.5: 


Form 1 the division into classes is very sharp, with Zone A dominating 
the costs. 

Form 2 less sharp division 

Form 3 no order of priority 
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Figure 5.5 The three forms of the ABC curve. 


5.3 RANK CORRELATION: SPEARMAN’S COEFFICIENT 
Ps 


Spearman’s coefficient is a measure of the relation between pairs of 
quantities that makes no assumption about their being normally dis- 
tributed and can be used with both continuous and discrete variables. 
All that is necessary is that the values of the variables in question can be 
put in order, either increasing or decreasing: this is the reason for the 
term ‘rank’. 


5.3.1 The method 


Let X, Y be the variables whose relationship we wish to investigate. 
The procedure is as follows. 


1. The values of X and Y are put in increasing order; suppose there are 
n, 

2. The first (smallest) value of X is given rank 1, the next largest rank 2 
and so on until all have been ranked; if there are equal values each is 
given the average rank for the group. This is repeated for Y. 

. For each pair (X;, Y;) the difference d; of the ranks is calculated. 

4. ps is calculated from 


eS) 


62d? 
(Siete Nigam 
ae ard! 


(a) Interpretation 
Ps has a value between —1 and +1; if 


Ps = 0 there is no correlation between X and Y 
= 1X and ¥ are strongly correlated and increase or decrease 
together 

—1 X and Y are strongly correlated and increase or decrease in 
opposite directions 
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A non-zero value can arise by chance even when X, Y are completely 
uncorrelated, and so a test for the significance of the value found is 
needed. If there are more than 10 pairs (X;, Y;) Student’s test is 
applicable, as follows. 

The quantity 


wig Ze 
= py? 
has a Student distribution with n—2 degrees of freedom; thus the 


hypothesis that pS is zero (i.e. that the variables are not correlated) is 
rejected atlevel aif T > t(n —2;-1-— a). 


5.3.2 Program for computing ps 
A program for computing pg is given in Fig. 5.6. 


10 CLS 

20 LOCATE 12,15:PRINT “RANK CORRELATION” 

20m LOCATE 13.1 PRINT “2. eee 

40 LOCATE 15,20:PRINT “SPEARMAN COEFFICIENT” 

50 FOR E=1 TO 3000:NEXT E 

60 PRINT:PRINT 

70 CLS:LOCATE 6,5:INPUT “NO. OF OBSERVATIONS _”,N 

80 LOCATE 9,5:INPUT “NAME OF FIRST VARIABLE _ :”,X$ 

90 PRINT 

100 LOCATE 12,5:INPUT “NAME OF SECOND VARIABLE Vs 

110 DIM RI1(N),R2(N),S(N) 

120 FOR I=1 TON 

130 CLS:LOCATE 4,5:PRINT “OBSERVATION 31; :” 

140 LOCATE 6,5:PRINT “RANK ACCORDING TO CLASSIFICA 
TION ”;X$;TAB(40) 

150 INPUT R1(I) 

160 LOCATE 8,5:PRINT “RANK ACCORDING TO CLASSIFICA 
TION ”;V$;TAB(40) 

170 INPUT R2(1) 

180 PRINT 

190 NEXT I 

200 D2=0 

210 FOR I=1 TON 

220 D2=D2+(R1(1)—R2(D) 2 

230 NEXT I 

240 R1=1—-(6*D2)/(N 3-N) 

250 PRINT:PRINT 

260 PRINT “OBSERVATIONS ”,X$,Y$ 

270 PRINT 


Figure 5.6 BASIC program for the calculation of the Spearman correlation. 
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280 FORI=1TON 

290 PRINT TAB(2);I,R1(1),R2(1) 

300 NEXT I 

310 PRINT:PRINT 

320 PRINT “SPEARMAN COEFFICIENT=__”: 
INT(10000 *R1+.5)/100;%” 

330 END 


Figure 5.6 (continued) 


5.4 ANALYSIS OF VARIANCE 


Correlation analysis takes into account only one factor at a time and 
therefore cannot measure the interaction of two factors. For this the 
method of analysis of variance is used, which enables us 


@ to study the simultaneous effects of several factors 
@ to reveal any interactions between different factors 
@ to optimize the number of observations needed 


It thus seems particularly well adapted to investigating the reasons why 
reject items are being produced. 


5.4.1 Mathematics of the method 
Analysis of variance is based on these assumptions: 


@ the factors affect only the means of the observations and not their 
variances; 

e the effects of the different factors are additive; 

e@ the residual variations (the ‘errors’) in the observations, after the 
effects of the factors have been taken into account, are distributed 
normally with zero mean. 


(a) The general linear model 

We consider the case of a quantity Y that is affected by two factors A 
and B which can interact; is Y; is a value observed for Y when A and 
B have values A; and B;, respectively, the model is 


V5 Se iy FB ey 


where a; and #; measure the effects of A, B respectively, Yi Measures 
the effect of the interaction between A and B, e; is the residual error 
and m is a constant. 

In order to study the interaction effect we must have repeated 
observations with the same values of A and B. The data can be set out 
as in Table 5.3 in which the following notation is used. 
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Table 5.3 Notation for two-factor analysis of variance 


Factor A Factor B 
B, By B; B, Total Mean 

a Yui Y 121 ijt Y 1st ies ee 

; Yin = Y120 Yijn Viege Sete alk 
A Deki Som y2jl Y2s1 Sree 

é Y2in Y22n Y2jn Y2sn vos Woe 
a Vit Yi2 Vij Vist 

: Yitn Yi2n Yijn Yisn i Ji 
A Yrit Yr21 Yr Versi 

j Yrin Yr2n Yrjn Yrsn bY ree Mito 
Total Va Y2. Va Ys. ae =. 
Mean Ya a Vis Vos a 
eS 


A,, Az, ... A, are the values of A 
B,, B,, ... B, are the values of B 


(These need not be numerical values. For example, they can label 
classes to which the factors belong.) 


Yi, (kK = 1, 2,..., n) is the result of the kth observation with A = 
Ai, BS B; 
Y; = yy is the sum of the values corresponding to (A;, B;) i.e. in 
i one cell of the table 


he 
= >'— is the mean of these values 
= 


yo = Se is the sum of the values corresponding to A,, i.e. the row 
ik sum 


= Vik . 

Y; = > is the corresponding mean 
pe 
hk 


i) Yj is the sum of the values corresponding to B,, i.e. the 
ok column sum 
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Y;.= Yeninr is the corresponding mean 

Y = Sy Y x is the sum of all the values 
i,j,k 

Y= Y__/nrs is the overall mean 


Yue Y_ is the deviation of an individual observation from the overall 
mean; the method considers the sum of the squares of all these 
deviations in order to attribute the total to the separate effects of A, B 
and the interaction AB respectively. It uses the identity 


Ge = ey = (Y ix xr Y;) ar (Y = Ya at (Y ;. <m YE) 
+(¥, Yet (Yj seen.) 


Squaring both sides and summing over i, j, k we get 
> (Yin iM alm ns> (Yi. = Yas a nr>, (Y zs aes 
ijk i i] 
a n> [(Y 5. = Yow) 3 (Y. ry Vale 
ij 


ots > Vin oa Yes 


ijk 
since all the sums of cross-products vanish. We can write this as 
SST = SSA + SSB + SSAB +-SSE 


where SST is the sum of the squares of all the deviations from the 
mean, SSA, SSB and SSAB are the contributions to SST that can be 
attributed to the effects of A, B and the interaction AB respectively 
and SSE is the residual (random) effect. This equation ‘analyses’ the 
total variance into the separate components; it can be shown that these 
components have a y* distribution with the following degrees of free- 
dom: 


SSA -—“w=r-1 

SSB ug = § — 1 
SSAB™=ti4zg = (r —1)(s — 1) 

SST uy = nrs — 1 

SSE pp = (n — 1)rs 


Dividing each of the sums of squares by its degrees of freedom gives an 
estimate of its variance; all this is summarized in Table 5.4. 

If A, B and the interaction AB have no effect on the observations — 
the ‘null hypothesis’ — the values SA? etc. are all estimates of the 
residual variance. It can be shown that in this case the ratios SA?/SE?, 


EY COLLEGE 


> 
a 
<= 
or 
ot 
AS 
ui <t 
— 
eae = 
c 
O 
ri 
(T i u) sd ( ) 
eS T — u)ss JOLIO [eNpIsoy 
ASS is? 
CN ee) Le ee 
=a ke 8)( feu BE ly ee wore 101 
avss g ( Xe Jima al A A qv I 
17" = as s © 1 SA) 
at = =e Ul = Ioj0e 
aoa I ce Am BA Ke ass g d 
1 Ge 2) : 
= -VS T-4 A- 'K US = 10}08 
vee AA A Ke vss V A 
4odda fO aouvidpa fo ajouysq 


wmopaad{ [0 saassaq 


sasonbs fo sung  jaaffa fo aainog 


OV[NUIIOJ SISeG :IOULIILA JO SISATEUL I0JDEJ-OM], p's BIGRI, 
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SB’/SE? and SAB?/SE? all have a Fisher distribution with degrees of 
freedom (v4, Vg), (Vg; VE), (Vag, Ve) respectively. 


(b) Significance test for the factors 

Given this last result and a table of the Fisher distribution — one is given 
in Appendix 4 - a significance test for the null hypothesis is easily 
constructed. 

However, if there is an interaction effect the separate effects of the 
factors are not additive and one of the hypotheses on which the method 
is based is not valid; in such a case we cannot measure the effects of the 
separate factors. It is therefore usual to start by investigating the 
interaction. 


(c) Testing the null hypothesis 
The Fisher function is denoted by F(v,, v2; 6). The significance test for 
the interaction is: if 


SAB?/SE? > F[(r — 1)(s — 1), rs(n — 1); 1 - a] 


then the hypothesis that the interaction AB has no effect is rejected at 
probability level w. The conclusion is therefore that there may be an 
interaction effect. 

Similarly, for the factor A, if 


SA?/SE? > F[(r — 1), rs(n -— 1); 1 - a] 
then the hypothesis that A has no effect is rejected at probability level 


aw, and so A may have a significant effect. Similarly for B, if SB’/SE* > 
Fis 1), (n=) Sa}: 


5.4.2 Program for two-factor analysis of variance 


A BASIC program is given in Fig. 5.7. When this is supplied with the 
experimental data as input it outputs the values of F(A), F(B), F(AB) 
for application of the Fisher test. 


Example 3 
The engineer responsible for the maintenance of a group of production 
machines wishes to know if the monthly cost of the substandard items 
produced by a machine is related to its age and/or its total annual 
production. 

Taking factor A to be the age of the machine and factor B its annual 
production the data are as in Table 5.5. The program gives these results: 


SAB?/SE? 2.9 no influence at the level w = 0.025 
SA?/SE? = 8.02 weak influence 
SB /SE? = 258.5 _ strong influence 


CLS 


REM “ANALYSIS OF VARIANCE - 2 FACTORS (R AND OC), 


DIFFERENT MODES, N REPETITIONS” 
LOCATE 10,8:PRINT “ANALYSIS OF VARIANCE 2 FACTORS” 
DIM Y(20,20,5) 

Y0O00=0 

YIJK2=0 

LOCATE 12,7:INPUT “‘r (number of modes of factor A)”;R 
LOCATE 14,7:INPUT “‘c (number of modes of factor B)’’;C 
LOCATE 16,7:INPUT “‘n (number of repetitions)”;N 
LOCATE 18,7:PRINT “Enter Y(i,j;k) one line at a time”’ 
FOR I=1 TOR 

FOR K=1 TON 

FOR J=1 TOC 

INPUT “x”’sx 

Y(,J,K)=X 

YIJK2=YIK2+Y(I,J,K) 2 

Y000= Y000+ Y(I,J,K) 

NEXT J 

NEXT K 

NEXT I 

YI=0: YI002=0 

FOR I=1 TOR 

FOR J=1 TOC 

FOR K=1 TON 

YI=YI+ Y(I,J,K) 

NEXT K 

NEXT J 

Y1002=Y1002+ YI 2 

YI=0 

NEXT I 

YIJ=0 

YJ=0: YJO02=0: YIJO2=0 

FOR J=1 TOC 

FOR I=1 TOR 

FOR K=1 TON 

YJ=YJ+Y(I,J,K) 

YIJ=YIJ+ Y(I,J,K) 

NEXT K 

YIJO2=YIJO2+ YIJ 2 

YIJ=0 

NEXT I 

YJ002= YJ002+ YJ-2 

YJ=0 

NEXT J 

SST=YIJK2—(Y000°2)/(N*R+*C) 


Figure 5.7 BASIC program for two-factor analysis of variance. 
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460 SSA=YI002/(C*N)—(Y000° 2)/(N *R *C) 

470 SSB=YJ002/(N*R)—(Y000 2)/(N*R *C) 

480 SSE=YIJK2—YIJ02/N 

490 SSAB=SST—SSA—SSB-SSE 

500 PRINT ‘Calculation of variances’’, 

510 SA2=SSA/(R-1) 

520 SB2=SSB/(C-1) 

530 SAB2=SSAB/((R-1)*(C-1)) 

540 SE2=SSE/(R*C*(N-1)) 

550 IF N=1 THEN GOTO 560 ELSE GOTO 590 

560 SE2=SAB2 

570 SAB2=0 

580 PRINT “NB If there are no repetitions, SE is estimated by SAB and 
SABn doesn’t exist” 

590 PRINT “SA2”;SA2 

600 PRINT ‘“SB2”;SB2 

610 PRINT “SAB2”;SAB2 

620 PRINT “SE2”;SE2 

630 FA=SA2/SE2 

640 FB=SB2/SE2 

650 FAB=SAB2/SE2 

660 PRINT “FA”;FA,“FB”;FB,“FAB”;FAB 


670 END 
Figure 5.7 (continued) 


Table 5.5 Data for Example 3 (p. 110) 


B1: 20000 B2: 50000 B3: 80000 


Als <3 yr Yun = 20 Y 121 = 35 Y 431 = 60 
Yin = 23 Yin oy) Y 132 = 58 
A2: 3-6 yt Yo = Y x1 — 36 Y 31 = 48 
Yo 19 Yoo =) Y 930 = 55 
AS: Ola. Vays 28 Y32, = 40 Y 337 = (70 
Y 312 she Y 329 = 45 Y 330 = 15 


5-5 EXPERIMENTAL DESIGNS OF TYPE 2” 


Designs of type 2” are designs in which there are n factors, each of 
which can be present at either of two levels. These levels need not be 
specified quantitatively, the only requirement being the possibility of 
distinguishing between the two - for example, for a lubricating oil, 
between high and low viscosity; for a transaction-processing system, 
between high message rate (say 3000 transactions per minute) and low 
message rate (say 500 transactions per minute). The aims of this design 
are as follows: 
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® to minimize the number of tests needed; 

to quantify the effect of each factor; 

®@ to quantify the residual variance, i.e. the ‘error’ observed after the 
effects of the known factors have been accounted for; 

@ to reveal any interactions between the factors. 


5.5.1 Designs without replication 


In a complete design, i.e. where every factor enters at both levels, there 
are 2” observations. Thus for three factors there are eight observations. 
Reduced layouts can be designed that require only 2""! observations, 
but for these the estimation of the effects of the various factors is more 
difficult. 


5.5.2 Designs with replication: the type 2”*” 


Here each experiment is repeated r times; this enables the effects of 
factors not taken into account explicitly to be measured, and so the 
residual variance can be estimated. 


Notation 
Upper case letters A, B, C etc. denote the various factors; the two 
levels at which each can be present are denoted by + (higher), — 
(lower). 

In Yates’s convention a lower case letter, e.g. a, indicates that that 
factor, e.g. A, is present at its upper level and that all the factors not 
mentioned are at their lower levels. This is shown in Table 5.6. It 
simplifies the specification of the various combinations. 


Table 5.6 Yates’ table of signs for 2° design 


Symbols Sign combinations for effects of factors and interactions 
A B Cc AB AC BC ABC 

(1) = = S + Ht re = 
a + - - - - + + 
b — + - - + - + 
ab + + - + _ - - 
c - - + - + - - 
ac ate = ats = aF = = 
be at rn an = = a a 
abc + a ae ay oe 5 ae 
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5.5.3 The 2° design 


(a) The Yates table 

In Table 5.6 the row labelled (1) gives the results of experiments in 
which all the factors are at their lowest levels (—). In row a, A is at its 
high level (+), B and C are low (—) etc. The table simply gives the 
rules for the signs for constructing the expressions of section 5.5.3(b) for 
the effects of A, B etc. 


(b) Estimates of the effects of the factors 
An estimate of the effect of any factor is the difference between the 
results of the observations made with that factor at its highest and 
lowest levels. Thus for the above design, 


effect of A = i[(a) + (ab) + (ac) + (abc) — (1) — (b) -— (c) — (bc)] 
and 
@ = 4(effect of A) 
effect of B = }[(b) + (ab) + (bc) + (abc) — (a) — (c) — (ac)] 
B = 4 (effect of B) 
etc. 


The effect of the interaction AB is obtained by the rule of signs in the 
table: 


effect of AB = 3[(1) — (a) — (b) + (ab) + (c) — (ac) — (bc) 
+ (abc)] 
aB = 3(effect of AB) 


(c) Treatment of the residual variance 

The residual variance represents the combined effects of all the factors 
that are unknown or otherwise not under control. If the investigation is 
to give meaningful results it is important that these residual influences 
are randomized as far as possible; otherwise there is a risk that some 
results will be unreliable and the residual variance will be increased. 


(1) Design with replications 
The estimate of the residual variance is 


g2 = Y (dijnu)? 
ijk = k i 1 


“aS 


where dix, = Yijtu — Y ij, and n—k—I1 is the number of degrees of 
freedom. 
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(ii) Design without replications 
The estimator is the higher order interaction; thus for the 2° design 


6 = apy 


(d) Three-factor model 
The model is 


Yin = YZ + @Xy+ Be; + aBX1;Xo; + Pax 
+ wY¥X 1; X3, + ByYXayjX3~ + ig 


where Y;, is the true value at the point i, j, k, i is the estimated 
value at that point and d;, is the error term. 


be ae re a 
diz Serves to estimate the variance about the regression line: 


(0, = =u 


where the denominator is 1 because with eight results there is only one 
degree of freedom for the variance. 
In fact, it is the second order interaction that provides an estimate for 


(a): 
(One = 8(dijx)? with ix = +tapy 


Ve. 
(6,)* = 8(taBy)’ 


(i) Significance tests b 
To test the significance of the coefficients @ 6, 7 and af we need 
estimates for their standard deviations 6,, 6g etc; in this case these are 

ve 
The test is based on the result that @/6,, B/6,, 9/6, and wB/Gp all 
have a Student distribution with one degree of freedom. 

By similar processes we can construct experimental designs of type 2” 
with n = 2, 4,5, 6,.... 


na 


Oy = Of = Gap 


= |x| = |aBy| 


Example 4 
The following results were obtained in a 2° experimental design. 
Ciuc b ab @ ac be abc 


£5 Vis 10 0 15 14 8 17 
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Table 5.7 is the application of the Yates table of signs to these values, 
with the sums corresponding to A, B, AB etc, together with the ratios 
of the estimated coefficients and standard deviations. The significances 
of these are tested by comparing them with the Student distribution with 
one degree of freedom. For a probability level 0.05 we find, from the 
table of Appendix 2, 1(0.05; 1) =6.3, and therefore we conclude that 
the only significant coefficients are those of A and AB. 
Thus the model for these results is 


w 


Y ix = 14.62 + 7.0X,; + 5.7X 1X2; 


Table 5.7 Yates’ table for Example 4 (p. 115) 


Combinations Results Total A B (Gy AB» A@— BEC “ABE 


(1) 15 ce = = = + “ + = 
a 18 a + — = - = a: + 
b 10 fe = =e = = + = + 
ab 20 “re “S = = + + = = 
c 15 ee = = 33 ap = — + 
ac 14 + a = “in - ap = = 
be 8 + = + + = _ + _ 
abc ‘Es a “5 oe ee a: 4° + oF 
= 117 PAR Pah AN 6) PS ae | 3 
Total Se = 1 Jove, “Asoo, 125" Oa 1oe 
Coef. 14.52. 2:62. 0.87, 1.12, 2.42, 0.62.70. 1250875 


5.6 GRAPHICAL METHOD: SCATTER DIAGRAM 


A plot of one variable against another gives a visual impression of the 
relation between them, and can give an idea of the degree of correla- 
tion. This is illustrated in Fig. 5.8. 


Example 5 

The following table gives the aptitude test scores for each of a number 
of employees and the production levels achieved by each over a certain 
period. The question is, is there any relation between these numbers? 


Employee 1 2 3 4 5 6 Z 8 9 
Test score 17 8 s I AS Ae tae 13 2 
Production 156:", 90°” 702-200" ALAN 100. S150) 150) 10 


i 
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Vf 
x 
x x * & 
ili oe No correlation 
X 
Vv VY 
a x 
. yi 
4 x 
Possible linear correlation 
Xx Xx 
Y Vf 
x 
x 
x ? ; % 
Possible non-linear correlation 
Xx Xx 


Figure 5.8 Examples of scatter diagrams. 


The graph of Fig. 5.9 suggests that the two are related. 


5 Opt Aes SPA), | OS 


Figure 5.9 Scatter diagram for Example 3 (p. 116). 


6 
Basic mathematics 


6.1 PROBABILITY: THEORY, DEFINITIONS 


The probability P(E) of an event E is a number lying between 0 and 1 
that measures, in some sense, the likelihood of the event occurring. It 
can be derived by enumerating all the possibilities: 


number of cases favourable to E 


E) = — oo 
ac total number of possible cases 


This gives the ‘true’ probability of E; the difficulty is to carry out the 
enumeration. The alternative is an experimental procedure based on 
observations: if a number of observations are made, in each of which E 
might occur, and if f,(E) is the fraction of observations in which E is 
observed to occur, then 


P(E) = lim f,(E) 


This is the ‘frequency’ definition of probability; the greater the number 
of observations, the closer the frequency approaches the ‘true’ proba- 
bility. 

As we said, 0< P(E) <1; conventionally, P(E) =0 corresponds to 
certainty that E will not occur, P(E) = 1 to certainty that it will occur. 


(i) Joint occurrence (intersection) of events 
The probability that both events A and B occur is 
P(A A B) = P(A)P(B\A) 


where P(B|A) means the probability that B occurs, given that A 
occurs. 
The general relation for n events A,, A,,..., A, is 


P(A, M A, (ales 2 M A,) a P(A,) P(A|A,) P(A3|A, A>) oon 
PAM A Arndt) 


If the events A, B are independent, meaning that the occurrence of 
either has no influence on the occurrence of the other, then 
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P(APRBY ="PCA)P(B) 


and in general, for n independent events, 


n 


P(A, G.AG Ole GA.) = 14s) 


i=l 
When the events under consideration are breakdowns of equipment the 
assumption of independence is usually valid. 


(ii) Occurrence of one or other of several possible events (union of 
events) 
P(A U B) = P(A) + P(B) — P(AF B) 


(as can be seen by considering the number of possible occurrences of A 
and B). 
If A, B are independent 


P(A U B) = P(A) + P(B) — P(A)P(B) 
and if A, B are incompatible (mutually exclusive) 
P(A 2M B) = P(A) + P(B) 


since then P(A M B) = 0. In general, for n events, 


OO aN ier SU ME = » P(A,) - py P(A;) P(A)) 
i i#j 
+ > PARP(A)PAD. 
i#j#k 
An important result that follows from the two-event case is as follows. 
1. If A, B are mutually exclusive but one or other must occur 
P(A) + P(B) = P(A U B) =1 


2. lf B is the converse of A, i.e. B is ‘A does not occur’, written 
B=not —A, or B= A, then 


P(AYS P(B) = 1 
and so 
P(=A) = 1 — P(A) 


6.1.1 Total probability, Bayes’ theorem 


Suppose a set of events E is made up of a number of non-overlapping 
Sais fs R=, Pen ie: 
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E={E} 


where E;™ E, = ¢ (the empty set) for all i, j with 14], and that an 
event B depends on at least one of the E;; then from what we have just 
shown 


P(B) = P(B\E;)P(E,) + P(BJE) + ... + P(BIE,) 


> P(BIE;) P(E;) 


Now if &; is any one of the E; 

P(E, OB) = P(E;|B)P(B) 
Therefore 
P(E; B) P(BIE)) 


P(E,|B) = P(B) = DE diraaae P(B) 


This is Bayes’ theorem: 
P(E;) P(BIE;) 


PEIB) = S PCEYPCBIED 


The importance of this theorem is that if the events E; are the possible 
causes of the effect B it enables us to calculate the probability that the 
cause was the particular event E; when B was observed to occur. 


6.2 PROBABILITY LAWS 


A number of laws describing probabilities are needed in discussing and 
measuring quality. They fall into two classes according to whether they 
concern discrete or continuous phenomena: the first relates to events 
that can be counted, such as the number of machine breakdowns during 
a given period; the second relates to measurements of physical quanti- 
ties such as length, weight or electrical resistance. 


6.2.1 Discrete laws 


(a) Binomial law 

The binomial law concerns sampling from a batch (of manufactured 
items for example) when the composition of the batch is not altered by 
the drawing of the samples; this will be the case when the sample size is 
small compared with the batch or when the sample is returned to the 
batch after examination—this is called non-exhaustive sampling. If 
samples of n items are drawn from a batch of N the condition for 
applicability of the binomial law is that n < N/10. 
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If p is the fraction of defective items in a batch (the percentage is 
100p) the binomial law B(n, p) states that the probability that a random 
sample of n will contain exactly k defectives is 


P(x = k)= 4 Pk = py" 


where 


i Es n! 

kk} k\(n — k)! 

The characteristics of B(n, p) are as follows. 
E(x) = np 


where E(x) is the mathematical expectation. This means that if the 
sampling is repeated many times the average number of defectives will 
be np. The variance is given by 


o°(x) = E(x — np)’ = np(1 — p) 


It is conventional to write q = 1-— p (q is the fraction of non-defective 
items) and so the variance is npq. 
The distribution is illustrated graphically in Fig. 6.1. 


n=10 p=10% 


inreZiico ty 4a eS 


Figure 6.1 Binomial law: probability distribution. 


The cumulative function F(k) is the probability P(x <k) that the 
sample will contain at most k defectives: 


P(x <k) => (") vi pyre 


This is illustrated in Fig.6.2. Tables are available for P(x =k) and 
P(x =k). 


Example 1 

A piece of electronic equipment requires four resistors; these are drawn 
from a large batch for which it is known that the fraction of defectives is 
5%. What is the probability that (1) three of the four will be defective 
and (2) at most three will be defective? 
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Figure 6.2 Binomial distribution: cumulative distribution. 


(1) The probability of three defectives is 
3) (0.05)3(1 — 0.05) = 0.0005 


(2) The probability of at most three defectives can be calculated by 
summing the probabilities for 0, 1, 2, 3 defectives, since these are 
mutually exclusive events; but since there must be either 0, 1, 2, 3 or 4 
in the sample the probability is P(x <3) = 1— P(x =4) = (1 —0.05)* = 
0.9999. 


(b) The hypergeometric law 

The hypergeometric law replaces the binomial law when the assumptions 
on which the latter are based cannot be made and the sampling does 
affect the composition of the batch. If, as before, the batch and sample 
sizes are N and n respectively, p is the fraction of defectives in the 
batch and g = 1 — p, the result is 


ro = i= (MP) (.N}/( 


The characteristics are as follows: 


E(x) = np 
(as for the binomial) and 
NG =H 
2 — 
OX) = a7 MPa 


The cumulative function 
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k 


F(k) = P(x S k)= 2 P(x = = /) 


=o 
This distribution has the same graphical appearance as the binomial. 


Example 2 
A batch of 25 items is known to contain five defectives. What is the 
probability that a sample of five will contain three defectives? 


| Cre Ae) () 2) es) 00307 


(c) The Poisson law 
This important law can be regarded in either of two ways: 


The answer is 


1. as a limiting form of the binomial when the batch size N is 
effectively infinite, the sample size n is very large and the probability 
p (here of defectives) is very small but the product np poe a finite 
value, m say; 

2. as a description of the occurrence events, e.g. if a machine might 
break down at any time in a certain interval and the average number 
of failures in an interval of that length is known, the law gives the 
probabilities of 0, 1, 2, ... failures occurring in the interval. 


The Poisson law is that if m is the average number of defectives in the 
sample of the size that is to be drawn (or the average number of failures 
in the interval of interest) then the probability that the sample will 
contain k defectives (or that there will be k failures in the interval) is 


P(x = k) = exp(—m)m*/k! 
The characteristics are 
E(x) =m 
and 
Ge (x) yi 
(the equality of the mean and variance is a strong characteristic of the 
Poisson distribution). 
Graphically this again is similar to the binomial. 


Standard tables of the Poisson distribution P(x =k) and of the 
cumulative distribution F(k) = P(x < k) are available. 
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Example 3 

It is planned to give a one-day demonstration of a certain machine for 
which the average number of breakdowns in a five-day week is known 
to be 10. What is the probability that it will not fail during the 
demonstration? 


Here the mean rate is m = 10/5 =2 breakdowns per day and so the 
probability or no failures in a day is 


P(x = 0) = exp(—2) 2°/0! = exp(—2) = 0.135 


6.2.2 Continuous laws 


Here we are dealing with variables measured on continuous scales and 
are concerned with such things as the probability that the value of a 
quantity x lies in a certain range or does not exceed a certain limit. 
Corresponding to the probability function P(x = k) of the discrete laws 
we now have p(x, x + dx), the probability that the value of a random 
variable X lies between x and x + dx, where usually dx is small. We 
write 


p(x, x + dx) = f(x) dx 


where f(x) is the probability density function. The corresponding 
cumulative function F(x), the probability that the. value of x does not 
exceed k, is 


F(k) = Px < k) =| f(x)ax 


where the lower limit for the integral depends on the range of values 
that x can take; this can be 0 or —© or some finite non-zero value. 


(a) The normal (Gaussian) law 
For a random variable with mean m and standard deviation o this is 


se 
fle) = syacenp|- S| 
with characteristics 
E(x) = 
and 
a(x) =o" 


Figure 6.3 gives a graph of f(x); it is symmetrical about x = m, i.e. 


f(m —x) = f(m +x). 
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mM 


Figure 6.3 Gaussian (normal) law: probability distribution. 


Quantities that are important in the use of this distribution are the 
fractions of the total population in intervals m+ ko, i.e. within k 
standard deviations on either side of the mean, for various values of k. 
The values are 


k=1 68.26% 
k=2 95.45% 
k=3 99.73% 


These are shown in Fig. 6.4. The cumulative function is 


F(x) =| f(@)d8 


where the lower limit is —~ because there is no restriction on the range 
of values that x can take. 


95.45% 
99.73% 


Figure 6.4 Gaussian (normal) law: critical ranges for distance from the mean. 
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It follows from the form of f(x) that F(+) = 1 (as it should be) and 
that F[(x — m)/o] =1— F[(m — x)/o]. F(x) is shown in Fig. 6.5. 

The integral for F(x) cannot be evaluated analytically, and so 
numerical methods such as Simpson’s rule have to be used. In practice 
standard computer programs are now used, or the standard tables that 
are available. The tables give the value as a function of the ‘reduced’ 
variable u=(x — m)/o, for which, for example, F(—u) =1-— F(u). 
Such a table is given in Appendix 1. 


Figure 6.5 Gaussian (normal) law: cumulative distribution. 


Example 4 

The resistances of a production batch of resistors are normally dis- 
tributed with mean 600m and standard deviation 120nQ. For a 
certain assembly the permissible upper and lower limits are 720 and 
420 mQ_ respectively. What percentage of the batch will meet this 
requirement? 


The fraction is 


720 
1 - (« = 600)* dx = F(u) — F(u;) 


ie 120V (27) we 20? 
where 
Uy = (720 — 600)/120 = 1 
u; = (420 — 600)/120 = —-1.5 
From the table we find F(1)=0.8413 and F(-1.5)= 


1— F(1.5) =0.0668 and so the fraction is 0.8413 — 0.0668 = 
0.7745 = 77%. 
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6.2.3 Central limit theorem, distribution of the mean 
This important theorem states that: 


if X,, X,, ..., X, are independent random variables with arbi- 
trary distribution laws then the distribution of the sum Y = 2, X; 
tends to the normal law as n increases. 


The justification of the use of the normal distribution in many applica- 
tions of statistical tests rests on this theorem. It is important that the 
influence of each X; is small and that all have more or less equal 
influences. 

Whatever the distributions of the X; 


aad-o-(Y) = 207 (X,): 

We are interested especially in the distribution of the mean of some 
property measured for each item in a sample of n. If these values are 
distributed normally with different means but all with the same standard 
deviation, the mean also is distributed normally, with parameters 


Polk Sis = o(x) 
n » pee) Vn 
A corollary is that the difference Z of two normally distributed 


random variables X, Y is distributed normally with parameters 


Y Aaa) i aide: dp Vaud LC el 99) 


6.2.4 The log normal law 


Here the logarithm of the random variable x is distributed normally; if 
as usual the parameters are m and o the probability distribution 
function is 


1 - (In x — m)? 
2Oe 


= x<0 


forx=0 


The characteristics are: 
E(x) = exp(m + 307) 
o?(x) = exp(2m + o7)[exp(o’) — 1] 


The cumulative function F(x) is computed by changing to the reduced 
variable u = (Inx — m)/o. Figure 6.6 illustrates f(x) and F(x). 
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Figure 6.6 Log normal law: (a) probability (b) cumulative distribution. 


6.2.5 The exponential law 


This is particularly relevant to the reliability of electronic equipment; 
like the Poisson law for discrete variables it describes events that can be 
considered as occurring at random, such as breakdowns. 

The probability density function is 


f(x) = exp (-Ax) 
for x = 0. The characteristics are 

Eoemas 

0° (x) = 1/7? 


The cumulative function is 


F(x) = |, f(@) ag = 1 — exp(-Ax) 


Figure 6.7 illustrates f(x) and F(x). 

This law is closely related to the Poisson law. Let A be the average 
rate at which certain events such as breakdowns occur, so that the 
average number in a period of length x is Ax; then if in the Poisson law 
the parameter m has the value Ax the probability of k breakdowns in 
the period x is 
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P(k) = exp (—Ax)(Ax)*/k! 
and the probability of no breakdowns in that period is 
P(O) = exp(—Ax) 


which is the exponential law. 


Figure 6.7 Exponential law: (a) probability (b) cumulative distribution. 


Example 5 

The average failure rate of a piece of electronic equipment is estimated 
as 1 failure per 100000h (10~° failures per hour). What is the 
probability that it will fail between 200 and 300 h in operation? 


Here the variable x of the law is the time t; the probability that the 
equipment will fail between 200 and 300 h is 


F(300) — F(200) = exp (—0.002) — exp (—0.003) = 0.001 


6.2.6 The Weibull law 


This is much used in reliability studies, particularly for mechanical 
systems. It has the advantage of being very flexible and can be adapted 
to the needs of a variety of circumstances. 

The probability density function involves three parameters: 


Zany \ $1 slereg As 
ee | 
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where x — y > 0, B is the shape parameter (a dimensionless number), 1 
is the scale parameter (dimension of the variable x (here, time)) and y 
is the location parameter, also of dimension x. The cumulative function 


is 
2 eae -(2=*) 
F(x) =1 exp| (== 


The characteristics are 


E(x) = y + nf + 1/8) 


where I denotes the gamma function (see Appendix 9) and 


o%(x) = rr ( 7 : nf E 3) | 


Figure 6.8 illustrates f(x) and F(X). 
6.3 CONFIDENCE INTERVAL FOR THE MEAN 


6.3.1 When the variance is known 


If X is the mean of a set of nm normal random variables X;, i = 1, 2, 
.., n, all of which have mean m and standard deviation o(X), i.e. 
2, X; 


n 


X= 


then X is distributed normally with mean m and standard deviation 
o(X) = 0o( X)/Vn. The symmetric confidence interval at probability level 
a is defined as 


X= uC = bo) — eo eek elm 4a) 


Vn 


where u(1 — 4a) is the value of the reduced normal variable for the 
argument 1 — 5a. 


6.3.2 When the variance is not known 


This has now to be estimated from the sample; the unbiased estimate is 
ee RE Meat 
OX ) = sonia tac 
ieee 


It can be shown that (X — m)/(o/Vn) is distributed as Student’s ¢ with 
n—1 degrees of freedom. It follows that the corresponding confidence 
interval is now 


Confidence interval for the mean 13] 


0 


Figure 6.8 Weibull law: examples of (a) probability and (b) cumulative distribu- 
tion. 


oO 


Vn 


<m<X + t(1— 4a; n -—1) — 


X — t(1 — fa; n - 1) V 
n 


6.3.3 Confidence interval for the mean time between failures 


This derivation is based on the assumption of truncated sampling with 
replacement of failed items. If n is the number of items tested, ¢ is the 
duration of the test and n is the number of failures recorded, then the 
estimate of MTBF is 
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MTBE = nt/r 


An interval AB, where A and B are the lower and upper bounds 
respectively, can be determined such that the probability that the MTBF 
lies in AB is 1— a— fB. A and B depend on a and f and are given by 


2nt 2nt 
A= same = ary 
X*(1 — B; 2r + 2) X*(a; 2r) 
For any particular case these bounds can be either calculated from 


tables of the y* distribution (e.g. Appendix 3) or read from the graphs 
of Fig’6-9. 


Example 6 
n = 1000, r= 4, t = 100h; a = 6 = 0.05. 


The estimated MTBF is 1000 x 100/4 = 25000 h. 
Lower limit = 0.39 x 25000 = 9750h 
Upper limit = 3.8 x 25000 = 95000 h 


6.4 LINEAR REGRESSION 


The problem here is that we have a set of pairs of observed values 
(x;, y;) — a cluster of points when plotted — and that we would like to 
represent them as well as possible by a straight line, i.e. to fit them to a 
model (Fig. 6.10) 


y=ax+b 


We have to find values fora and b. The procedure used in regression 
analysis is the method of least squares. This finds the values of a and b 
that minimize the sum of the squares of the differences between the 
observed values y; and the values of ax; + b given by the model. These 
differences y;— ax;—b are the ‘errors’ e; in the model’s predictions. 
The method minimizes E = >’ (e;)*. 

We have : 


a= 0% — ax; — b)? 
= > (y)? — 2a Dd xy; — 2nby + a? DS (x)? + 2nabz% + nb? 


where v is the number of pairs of (x;, y;). The values of a and b that 
minimize E are given by the solutions of 
Oe: ie OE 
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Figure 6.9 (a) Confidence intervals, truncated sampling: factors by which the 
estimated MTBF should be multiplied to give upper and lower limits at a stated 
level of confidence. 
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Figure 6.10 Linear regression. 


Differentiating and solving the resulting linear equations gives 
DLV, eV 
Pion ee Ea e, 
= (x;)* — nx? 
This gives the ‘best’ line in the least-squares sense. 

The method can be used for other models. Thus for the quadratic 
model y = ax* + bx +c we find E, the sum of the squares of the errors 
e;, as before and solve JE/da =0, JE/db =0 and JE/dc = 0 for a, b 
and c. 

Standard programs are available for fitting straight lines and other 
curves to statistical data. 


b = Ver ax 


Example 7 
For the values 


420" 40. 110" 60. 190 
y 60 90 120 140 170 
the method gives y = 0.62x + 48. 


Exercised 


(Solutions are given on pp. 148-153) 


i 


An automatic machine makes spacer bars whose length must be 
between 37.45 and 37.55 mm; the lengths produced have a normal 
distribution with mean 37.50 mm. 


(i) What must be the standard deviation if 998 out of every 1000 
bars are to be acceptable? 

(ii) A random sample is drawn from the production and the lengths 
are measured. What must be the size of this sample if the mean 
of the lengths is to lie between 37.495 and 37.505 with proba- 
bility 0.95? 


You are given that if z is the reduced central normal variable 
POO Sz = 1.96) = 0.475" 
P(O S z < 2.05) = 0.480 
P(O S z = 3.10) = 0.499 


. An automatic machine makes items whose weight is distributed 


normally with mean 0.90 g and standard deviation 0.06 g. 


(i) What is the probability that the weight of an item chosen at 
random lies between 0.84 and 0.99 g? 

(ii) How many items of weight less than 0.81 g can one expect 
there to be in a batch of 5000? 

(iii) The items are packed into boxes, 100 to a box, by another 
machine; a certain number of boxes are chosen at random and 
the mean weight w of an item is found for each box. 


(a) What are the mean and standard deviation of w? 

(b) What is the probability of a measurement deviating by 1% 
from this mean value? 

(c) What are the 95% confidence limits for the mean item 
weight in a box of 100? 
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(d)The mean weight of the items in a box of 100 chosen at 
random is found to be 0.88 g. Can the box be considered as 
representative at the 1% level? 


3. The lifetime, in hours, of an electric light bulb is a normal variable 
of mean M and standard deviation 20. A test of a sample of 16 
gives a mean life of 3000h. Find the 90% confidence interval for 
M. 


4. We wish to calculate the probability that a structural element will 
break under deflection. The method employed is what is called the 
R/C - resistance/constraint — method and involves constructing the 
random variable R — C, from which the probability can be deduced. 

The resistance R is normally distributed with mean R = 28 and 
standard deviation Og = 2. 
The constraint C is normally distributed with mean C = 25 and 


standard deviation 0¢ = 1.5. 


(i) Give the theoretical calculation. 
(ii) Find the probability of breaking. Units are daN/mm? (deca- 
newtons per square millimetre). 


5. A machine makes items whose diameter X is a normal random 
variable with mean 32 mm and standard deviation | mm. 


(i) What is the probability that the diameter of an item is less than 
30.5 mm? 
(ii) What is the probability of a diameter between 31 and 33 mm? 
(iii) As a control, samples of 20 are taken at regular intervals and 
measured; if the mean diameter found is X 


(a) what is the probability distribution of X? 

(b)in what interval [a,b] must X lie if the machine can be 
considered to be correctly adjusted with a probability of 
0.99? 


6. A factory makes a certain item in large numbers. This is done in 
two stages, in the first of which a defect A can appear and in the 
second a defect B. Experience has shown that 2% of the items show 
defect A and 8% show B. 


(i) Find the probability that an item chosen at random 


(a) has both defects, 

(b) has at least one of the defects, 
(c) has one and only one defect, 
(d) has no defect. 
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(ii) A sample of 200 items is taken and the number X showing 
fault A is noted. 


(a) X is regarded as a random variable with Poisson distribu- 
tion. What justification is there for this? What is the 
parameter? 

(b) What is the probability that 10 items in the sample of 200 
will show fault A? 


(iii) A sample of 300 is taken and the number Y showing fault B is 
noted. 


(a) Y is regarded as a normally distributed random variable. 
What are the parameters? 
(b) Calculate prob(Y < 24) 
prob(20 < Y < 35) 
prob(Y < 30, given that Y > 24) 


7. The objective of this problem is the investigation of the performance 
of a lathe, one of the machines in a company’s mechanical work- 
shop. The problem is in two independent parts. 

The lathe turns shafts to a nominal diameter of 24 mm; the actual 
diameter is a normal random variable of mean 24 and standard 
deviation 0.02 mm, and the tolerance limits are 23.95, 24.05. 


(i) How many good items will there be in a sample of 1000? 

(ii) As a control, 20 items are chosen at random and the mean 
diameter X is found; what is the interval [a, b] in which this 
must lie if the lathe can be regarded as correctly adjusted with 
probability 0.99? 

(iii) A check on 20 items gives the following values: 


Diameter Number Diameter Number 
23.93-23.95 1 24.01-—24.03 8 
23.95—23.97 41 24.03-24.05 Z 
23.97-23.99 1 24.05—24.07 0 
23.99-24.01 i! 


(a) Find the mean and standard deviation of this set of values. 
(b) What conclusion do you draw from these measurements? 


8. A quality assurance service decides to impose controls on dimen- 
sions X and Y of a product, using control charts for the mean and 
standard deviation. It has to estimate the population standard 
deviation o, and has the following observations: 


15.050 
15.055 
15.056 
15.060 
15.062 


15.049 
15.058 
15.061 
15.046 
15.054 


15-052 152056 
15.061 15.058 
15.070 15.049 
15.061 13:059 
15.07 15.080 


30.021 
30.025 
30.031 
30.026 
30.027 


(i) Give the control limits for the two machines. 


30.035 
30.032 
30.035 
30.032 
30.035 


30.020 30.021 
30.019 30.025 
30.022 30.024 
30.018 30.023 
30.017 30.040 
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(ii) Set up the required control charts, for a sample size 5. 


A study of the way the dimension X changes with time enables the 
effect of wear to be estimated, and consequently the time the 
machine tool can be allowed to run before there is a risk of rejects 


being produced. Observation gave the following: 


Operating time AA 


1 min 
8 min 
35 min 
50 min 
70 min 
80 min 


15.000 
15.010 
15.040 
15.070 
15.085 
15.090 


By fitting a straight line X =at+b to these values find the 


intervals at which the machines should be re-set. 


. A manufacturer of food products is looking for the most economical 
way of putting a powder into cartons. Trials using an old machine 
gave the following results: 
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Sample no. Weight (g) 
1 510 522 520 514 516 
Z 514 516 512 514 520 
3 516 514 518 512 514 
4 510 508 520 516 514 
5 516 504 512 516 22 
6 518 510 512 518 514 
i SZ 512 508 512 520 
8 514 518 514 512 516 
9 518 S17 515 514 510 
10 520 515 514 508 513 
11 518 514 516 312 516 
ily. 510 511 Say 512 510 
13 520 518 506 518 510 
14 518 bills) 516 2 512 
15 ya 506 510 522 522 
16 516 516 514 510 516 
17 516 510 516 520 522 
18 514 512 514 518 S12 
19 524 502 516 520 508 


20 514 514 524 516 518 


Here each observation is the total weight of the box and its filling; it 
is known that the box weight is distributed normally with mean 64 g 
and standard deviation 0.05 g. 


(i) (a) Can such measurements be used to keep the production 
under statistical control? Use either a goodness-of-fit test 
or Henry’s line. What are the relative advantages and 
disadvantages of the two methods? 

(b) Give the different estimators that can be used for o. 


(ii) The customer requires that the average content of a box shall 
be not less than 445 g and that at most 2% can be more than 
7 g short. 


(a) What percentage of boxes will contain less than the mi- 
nimum weight stipulated by the customer? 

(b) If the mean is changed what must its new value be to stay 
within the terms of the contract? 
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(ili) Set up the control chart that must be used for ensuring that 
(ii) is satisfied, assuming a sample size of 5. 
Samples of the output were taken at different times and the 
following weights were found: 


at 10, 00un S162, cR514)., AOL 2 tet 20s ly POLS 
10:30 Wbl4A wots nln 520, wag od 2 
11200 44 S10cj oa St6e \AS20 fc 5 1A aa 22 
11,530 55225. 90S 2 pee 51 2d eo D020 le OS 


Incorporate these in your control chart. 

(iv) What would be the cost to the manufacturer of adhering 
strictly to the customer’s conditions ((1i)(a) above)? The 
annual sales are 5 hundred thousand boxes of 445 g, priced at 
£1.50 per kilogram gross weight. 

(v) In response to invitations to tender for new packaging ma- 
chinery the manufacturer has these proposals: 


Machine A: o, = 1 g guaranteed 

cost £100 000 

estimated life 10 yr 

residual value after 10 yr £9500 
Machine B: og = 3 g guaranteed 

cost £55 000 

estimated life 10 yr 

residual value after 10 yr £5000 


The current value of the existing machine is zero. 
What should be the decision, on economic grounds? 


10. A workshop has 30 workstations, each equipped with a machine A 
and a machine B. The two machines operate completely indepen- 
dently of each other. The probabilities of breakdown in the course 
of a working day are P, = 0.1, Pg = 0.03. 

(i) For an arbitrarily chosen workstation during a working day 
find: 


(a) the probability P, that both A and B break down; 

(b) the probability P, that at least one machine breaks down; 

(c) the probability P; that A, and A alone, breaks down; 

(d) the probability P, that one, and only one, machine breaks 
down. 


142 Exercises 


(ii) For an arbitrarily chosen group of 10 machines of type A and 


(iii) 


(iv 


(v 


) 


—_— 


a single working day let X be the number that break down 
during that day. Give the probability distribution of X. 


(a) Find the probability that exactly three of these 10 ma- 
chines fail. 

(b) Find the probability that at least one fails. 

Give these probabilities to an accuracy of 1074. 

Let Y be the number of workstations at which one and only 

one machine breaks down during the working day. Justifying 

the use of the Poisson law, calculate: 

(a) the probability that four of the 30 stations will experience 
a breakdown of one and only one machine; 

(b) the probability that at most two stations will experience a 
breakdown of one and only one machine. 


Give these probabilities to an accuracy of 10~*; for (b) take 
P,=0.12. 

Two dimensions x, y are to be controlled for an item 
produced by the workstations, and the limits are x =6 + 0.01, 
y =3+0.01. Measurements on the production of a station 
give the mean and standard deviations: E(x) = 6.003, 
0,= 0.005; E(y) = 3.002, o, = 0.005. 

Assuming the x, y are normally and independently dis- 

tributed, what percentage of rejects can be expected in a day’s 
output? 
Because of the high cost of maintenance of machine A the 
management decides to reconsider the maintenance policy. 
From records of the numbers of fault-free working days over a 
year it appears that the lifetimes of these machines follow a 
Weibull law with parameters y = 80, n = 110, B =2.2. 

Let the random variable t be the number of fault-free 
working days before the first failure. On Weibull paper show 
the cumulative distribution function F for this variable, for 
t € [100, 250]. 


(a) Find the mean number of fault-free working days. 

(b) Hence find the probability of not failing during a period of 
this length. 

(c) Using the graph of F find the number of days after which 
30% of the machines will have had their first failure. 


11. A search through the maintenance file has revealed 27 failures of 
pumps of a certain make during the past 2 months. The times of 
fault-free operation, in increasing order, are as follows: 


i. 


13. 
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26 132 245 h 
43.4 145 247 
58.8 151 282 
68 171 295 
80.5 180 307 
94 192 320 
98 202 350 
112 210 335 
125 220 474 


(i) Using Weibull paper and the table of the Weibull law provided 
(Appendix 5) find the MTBF: use the method of mean ranks to 
determine the cumulative relative frequenciees F, = (2;n;)/ 
(N+ 1): 

(ii) Hence determine the maintenance routine best suited to this 
product. Justify your choice. 


During the period 1.9.89 to 1.9.90 a certain machine has been in 
operation for 1205 days and has failed 10 times. The intervals 
between successive failures were as follows (in days): 


AA NOS S2". 39-108 299° 57° 255. 151 49° — 


(i) Tabulate these values in a manner that shows the numerical 
order of the failure, the times between successive failures in 
increasing order and the accumulated good time. 

(ii) Let R(t) denote the probability that the machine will run for a 
time ¢ without failing. Using log-linear paper, justify the 
representation of R(t) by the exponential law. 

(iii) Find the MTBF for this machine. 


In order to obtain information concerning the reliability of a certain 
machine the performance of the example installed by a customer has 
been studied, and the following intervals between successive failures 
have been recorded: 


1704225, 260 300:°132007330* 390” :490 h 


(i) Find the Weibull model. 
(ii) Find the maintenance interval corresponding to a 50% risk of 
failure. 
(iii) (a) Draw the curve of 


= opiiqil 
wb! 


for t = [200, 400]. 
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14. (i) 


(il) 


(iii) 


(iv) 


(b) Interpret your results. 
The following lifetimes were recorded for a certain device: 
1860 2500 2900 3600 3950 5100 6300 h 


(a) Are these described by a Weibull law? 
(b) If so, what are the parameters? 

(c) Compute E(f). 

A control scheme is defined as follows: 


n = 20, A =2, R = 3, AQL = 2%, P, = 10% 


a, B denote the risks borne by the supplier and the customer 
respectively. 


(a) Compute P,, a, f. 
(b) Draw the efficiency curve for the scheme. 


In the system of the diagram element 1 is working; if it fails it 
is switched out by S and element 2 is switched in. The 
following assumptions are made: 


@ switch S§ has reliability = 1; 

@ the repair time follows an exponential law with parameter 
N; 

the time of fault-free operation of each element follows an 
exponential law with parameter A; 

e att =0, P(0)=1, P(1) =P(2) =0. 

We are interested in the availability of this system. 


(a) Draw the Markov chain. 

(b) Give the state-change equations. 

(c) Apply the Laplace transform and give the resulting equ- 
ations in matrix form. 

(d) Compute the availability. 

(e) What is the limiting value of this as t > ? 


The question concerns an aircraft blind-landing system. 
The risk of failure for any single flight must be less than 
10-°. An aircraft has 10 essential (electronic) systems, for 
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each of which the risk of failure in any one flight must be less 
than 107’ (this is the risk of the complete blind-landing system 
failing). The landing system is in operation only during the 
6-min approach to landing. A total of 1000 components are 
involved, for which the average value of A is 10~°/h. 


(a) Find the probability that the system will fail during a flight. 
(b) Does the system attain the specified risk level? 
(c) If not, what should be done to achieve this? 


15. To investigate the reliability of certain machines, a number of which 
are installed in a workshop and operated under the same conditions, 
the number of days of operation before the first failure have been 
recorded for 10 of the machines: 


SO 1106886 100 OE 1202 9413554. 

(i) Using Weibull paper, justify the representation of these times 
by a Weibull law with y = 50. Find the other parameters f, n 
and the MTBF. 

(ii) Find the running time for which the risk of failure is 50%. 

(iii) Given that the failure rate is A(t) = (B/n)[(t — y)/n]*"', evalu- 
ate A(t) for t= 60, 80, 100, 120, 130 and plot be values. 
Interpret the result. 


16. In a study to find estimates F(t), R(t) and A(t) for the reliability 
functions for a certain device the times to failure of 55 examples 
were recorded, with the following results: 


Time interval No. failing in 
(hy this interval 


U 


0-500 
500-1000 
1000-1500 
1500-2000 
2000-2500 
2500-3000 
3000-3500 


a 
SNIOnNN OO CO W 


(i) Give the estimates. 
(ii) Draw the graph of A(t). What do you conclude from this? 


17. In a food products factory the conditioning operations and the 
handling of the products preparatory to despatch are automated. 
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This automation requires the inclusion of 80 extra devices, all of the 
same type and all subject to very heavy use — 1.5 x 10° handling 
cycles each per month. The quality engineer, in order to set criteria 
for reliability, has recorded the failures of a sample of 50 of these 
devices, with the following results: 


Handling No. of 
cycles failure 


=) 


OPS TS 
1:5-=3 
3 - 4.5 
4.5- 6 
6 =-7.5 
7.5- 9 
9 -10.5 
10.5-12.0 
12.0—13.5 
13.5-15 
15 -16.5 
16.5-18 
18 -19.5 
19.5-21 
2 S755 
23.5-24 
24 =25.5 
25.0=2/ 


POrRrPFNNWANKANADNANHWN PH 


In general, the Weibull law is applicable to this type of equip- 
ment. 
(i) Show how the special Weibull paper can be used to test the 
validity of the Weibull law in this case. Expressing the law in 
the standard form 


noe co (aha 


find the values of 6, y, 7. 

(ii) Having found these values you wish to be assured that the 
Weibull law really does describe the data, and for this you apply 
the x’ test. How do you go about this? Give the full details of 
the calculation, and apply the test at the level wa= 0.05. [Table 
of y’ provided in Appendix 3] 
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18. A company that manufactures shipboard electric motors for the 
Navy must, without fail, ensure that the noise level of its products 
does not exceed a certain value (given in decibels), so that it is not 
detectable by enemy submarines. Investigation by the research 
department has shown how to deal with the well-known causes of 
excessive vibration (e.g. unbalanced rotor), but there are other 
factors whose effects are less easily identified: 


Factor A: viscosity of lubricant, Al = low, A2 = high 
Factor B: method of tightening the cover, 

B1 = uncontrolled, B2 = with a torque wrench 
Factor C: play in bearings, C1 <6 um, C2 >6 wm 


To investigate these a 2 X 2 X 2 experiment was undertaken. The 
Yates table for the results was as follows: 


(1). £70 
ee 145 

b 6) 
ab 167 
CrAs9 
ac 145 
be 78 
abc 167 


Use the Student test with confidence level 90% to find which factors 
have a significant effect. 


19. (i) Find the overall reliability of the following system, giving 
R(i) = 0.9. Compare the result with that for a system without 
redundancy. 


(ii) Repeat the calculation for the following system, with R(i) = 0.9 
also. 


Solutions 


1 

(i) = 0.016 

(ii) n = 40 

2 

(i) 77.4% 

(ii) 334 = 

(iii) (a) X , = 0.90 g; o, = 0.006 g 
(b) P = 0.88 
(c) No 

3 


—1.64< M <1.64 


4 

Gi) _P rupttre="PUR = 6) PRK —C'=— 0) po gen 
R-—C is normally distributed with mean m= R-— C, standard 
deviation 
o = (0% + o%)” 


whence 


(PIR = <0) = Pi <= —S) 


(ii) P(U < 3/2.5) = P(U < -1.2) 
P(R-C <0) =0.12 


5 
(i) 6.68% 
(ii) 68.26% 
(iii) (a) |X is normally distributed: 


6 


_ Oy 
Nixes 
ea 
(b) The interval is AB 
Sl AQ oN 3257 


(i) (a) P = 0.16% 


(b) P=9.84% 
(c) P=9.68% 
(d) P=9.52% 


(ii) m =4 


P(A) = 0.53% 


(iii) 7 is large 


“4 


m = np = 24 

o = 4.7 

P(Y < 24) = 0.41% 
P(20.< Y < 35) = 264% 
P(25 < Y < 29) = 27.22% 


(i) 986 


(ii) [23.948 — 240.56] 
(iii) (a) X = 24.006; o = 0.028 


(b) Adjustment is correct 


8 
(i) For X: 
oO =16.4 KO 
Li, = 14.941 
Lis = 14.944 
Lig = 14.989 
L., = 14.986 
Limits of warning zone: 
Tl — 66, =31.4'x.10" 
Ge se cls 
Time between re-setting: 
oS 
ime 27 min 
9 
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(i) (a) 7° is a more rigorous test for goodness-of-fit to the normal 


law; however, m and o have to be estimated in order to 
apply this, and here Xx» = 450 g, o, = 4.97 g, Wee xX» +X, 


The results are Petcsion 
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Ls 1 Ae 
(b) Ody ee poy Us Xx) 
j=m 
and d= — > o,-, 
j=l 
Or = (0%, + 03)” 
ss 438 — 445 
(ii) (a) U,= Gases 1.4 
Thus 8% are outside the tolerance limits. 
(b)) X= 4480 
(iii) Control chart for the mean (product + box) 
Ly, = 525.6 
Lo = 522.7 
Li, = 501.3 
Lic = 498.4 
Control chart for product + box 
£6. =e LG 
Ly, = 7.46 


X12. 28 2o-— 35 
Xe = 514, Si nO, Actes 
X;=511.4g 6,=47g 
X4 = 509 g 6, = 7.92 g 
Take another sample immediately after the fourth, to check the 
stability of the standard deviation. 
(iv) Base the calculation on 8% of boxes being outside the tolerance 
limits: number of boxes rejected = 5 x 10° x 8% = 4 x 104 
: 4 ae 
Cost: 4 x 10° x 445 i000 £267 000 
(v) Machine A, loss 0% 
Machine B, loss 1%, cost £33395. 
Choose machine A. 


10 
(i) (a) 0.003 
(b), 0.127 
(c) 0.097 
(d) 0.124 
(ii) (a) 2.96% 
(b) 74.8% 
(iii) (a) 0.19 
(b) 0.87 


(iv) 14.3% 
(v) (a) 179.4 days 
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(b) 0.46 
(c) 150 days 


‘| 

May = 0: Pis16; n= 225 
B=1.6=> A = 0.8966 
MTBF = 0.8966 225 = 201h 

(ii) 6 > 1 > preventive maintenance 


12 
Graphical solution 


(ii) 
(iii) MTBF = 120 days 


13 
(i) R(t) = exp[—(1/370)] 
(ii) ¢ = 330 


(iii) Rate increasing in the ‘aged’ period 


14 
(i) (a)  Weibull’s law with y > 0 
(b) = 2900h; B =1.5; y = 1200h 
(c) E(t) = MTBF = 1200 + 2597 = 3797h 
(ii) (a) P,; =NQA 


B = 0.676 
(b) Graph as shown 


100% 
99.28 


67.60 
Pa 


0) 2% 10% P% 


@) eee) 
Ww 


uw 


(iii) (a) 
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(b) Pot) = APo(t) + uP, (t) 
Pi(t) = APo(t) — (u + A)Pi(t) + MP Q(t) 
P(t) = AP;(t) — wP2(t) 


(s + A) (—w) 0 0 1 
(0) | A (stat+p)  (-w) | [1] = [0 
0 (—A) (s + py} -\2 0 


[s is used here for the parameter of the Laplace transform; in the text 
and the table, p] 


2 
ue + Au 
d A(t) = ————-_—_ 
(d) (t) eect 
* 2 Cyréexp(4r2t) '— see 
M+wt+hu Cre @, 
244 
(e) u : uu 
V+uw + du 


(iv) (a) 0.9990 
(b) = ADSM’ 
(c) Active redundancy 


its) 
(i) y = 50; B = 1.5; n = 32 
MTBF 95 h 
(ii) tsq% = 88 h 
(iii) t 60 80 100 120 130 
A(t) O:012" ©0023" ©0:032- 07039" "01042 
i (t) 
t 
16 
Time No. of No. of 
interval failures survivors - - 
i Len) Ue Uae) to ¢; R(t) F(t;) A(t;) 
0-500 8 35 1 0.054 1.09 x 10-4 
500-1000 8 Sy) 0.94 0.14 3.07 > 0m: 
1000-1500 10 44 0.8 0.18 Aes Se = 


—— eS 


Time No. of No. of 
interval failures SurviVOrs : 
bi — G41 bee bee to f; R(ti) 
1500-2000 12 34 0.618 
2000-2500 7 22 0.4 
2500-3000 8 ‘LS 0.27 
3000-3500 v7] 7 O27 
i (t) 
0 
t 
Increasing 
Adopt preventive maintenance 
‘Ageing’ period 
17 
(i) (a) v= On = 1338 X 10° eycles; B = 2.3 
(b) Test 7 
Number of classes ~ 8; E = 0.47 
8-3-1005 = 9.49: accept the model 
18 
@= 155 
B : = 19.3 
ap = 30.3 


with Y.. = 140.5 


19 

(i) R(s) with redundancy = 0.88 
R(s) without redundancy = 0.66 

Gi) R(s) = 0.96 


F(t) 


0.218 
0127 
0.145 
0.127 


Solutions 153 


A(t) 


To elO4 
636 <10-* 
10.4 x 10-4 
20 4x 1054 


Appendices (Tables) 


1. GAUSSIAN (NORMAL) DISTRIBUTION 


The table gives the cumulative distribution function F(u) = Jf“. f(v) dv 
where f(v) is the Gauss or normal law (p. 125). 

If x is distributed normally with mean wu and standard deviation o 
then u is the reduced normal variable, u = (x —)/o. F(u) is the 
probability of finding a value less than or equal to wu. 

The table is for positive values of u only; values for negative u are 
found from 


F(—u) = 1 — Flu) 
eg F(—0.94) = 1 — F(0.94) = 1 — 0.82639 = 0.17361 


The second table (p. 156) gives 1 — F(u) for large values of u. 


0L9L6'0 
79 0L6'0 
L7£96'0 
67 7560 
80776 0 
68 1£6'0 
pL Li60 
Ly 1060 
86 7880 
v1 7980 


16 8¢8°0 
LZE18 0 
y7S8L'0 
06 7SL'0 
Or ZCL'0 
€6.L89'°0 
eL1s9'0 
607190 
ce SLs'0 
98 s¢S'0 


ST 9L6°0 
56 696'0 
97 796'0 
7S €$6'0 
66 76 0 
9S 0¢6'0 
179160 
CL 6680 
00 188°0 
£6 6$8'0 


9 9¢8°0 
LSOT8°0 
O¢ 78L'0 
SLISL'0 
70 61L'0 
6€ 789 '0 
€0 8790 
9Z019'0 
7 LLS'0 
gg ics 0 


BS SL6'0 
97 6960 
79 196'0 
vs 756 '0 
6L 1760 
77 676 '0 
997160 
96 L680 
00 6L8'0 
69 LS8'0 


96 €£8°0 
58 L080 
9€ 6LL'0 
LS 8bL'0 
99 SIL'0 
78. 089'0 
Te 7r9°0 
Zr 909°0 
os L9S°0 
06 L7S°0 


00 SL6'0 
9¢ 8960 
08 096°0 
7S 1860 
79 0V6'0 
98 L760 
60 €16'0 
L1 968'0 
86 9L8°0 
cy SS8'0 


Ly 1¢8'0 
11 S08'0 
LEQLL'0 
Le spl 0 
977IL'0 
v7 LLO'0 
8s 0r9'0 
LS 7090 
9¢ £950 
76 £7S 0 


Ir vL6'0 
78 L960 
76 6S6'0 
€S 0S6'0 
ev 6£6°0 
Ly 960 
6r 1160 
SE 7680 
€6vL3°0 
v1 ¢s8'0 


76 878 0 
ve 7080 
LEeLL 0 
ST tL 0 
78 80L'0 
$9 £L9'0 
€8 9¢9'0 
IL 86S 0 
79 6SS ‘0 
7661S'0 


18 ¢£L6 0 
71 L960 
L06S6'0 
0S 6760 
TT 8LO 0 
LO S76'0 
88 6060 
TS 7680 
98 7L3'°0 
€8 0S8'°0 


6€ 9780 
6S 66L'0 
SE OLL'0 
168¢L 0 
Or SOL'0 
€00L9'0 
LO¢ee9'0 
78 76S 0 
L9SSS'0 
colts 0 


07 €L6'0 
ge 996'°0 
61 86 0 
Sv 876 0 
66 9¢6'0 
79 £760 
77 806 0 
69 0680 
9L0L8'°0 
0S 8780 


18 €78'0 
€L96L'0 
Te LOL'0 
$9 SeL'0 
76 1OL'0 
or 999°0 
0¢ 6790 
66 06S 0 
ZLISS'0 
L611S'0 


LSTLO'0 
79 596 0 
87 LS6'0 
ge Lv6'0 
yLSc6'0 
07 7760 
8S 9060 
LL888°0 
79 898°0 
v1 980 


1Z 178°0 
68 £6L'0 
v7V9L 0 
Le7eLl 0 
Ly 369'0 
9L799'°0 
7S S790 
90 L8S'0 
9LLvS'0 
86 LOS 0 


700 


€6 1L6'0 
68 796'0 
LE 9S6'0 
0 976'0 
sr 7c6'0 
€L 0760 
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2. STUDENT t DISTRIBUTION 


1 hieaiy 


Values (positive or negative) of ¢ having probability a of being 
exceeded 


F(t) = 


42\~vt1/2 t 
hee 4) F(t) = I F(t) dt 
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3. x* DISTRIBUTION 

XEN orks ms X,— m> 

Va ee 
O71 0» Gy 


ee 1 , v/2-1 ty 
FOO) = sara) 2)" exp| e| 


For v>30 the quantity (2y’)/? —(2v—1)'” can be taken to be a 
reduced normal variable. 
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4. THE F (FISHER-SNEDECOR) DISTRIBUTION 
= xi/ V1 
13/ V2 
Xi. 16 have v,, vy degrees of freedom respectively. The table gives values 
of F having probability a of being exceeded. 
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5. MEAN TIME BETWEEN FAILURES FOR A SYSTEM 
FOLLOWING THE WEIBULL LAW 


The table gives the values of A, B where the mean and standard 
deviation of the MTBF are found from 


mean = An + y standard deviation = Bn 


B A B 6B A B 


1.80 0.8893 0.511 

1.85 0.8882 0.498 

1.90 0.8874 0.486 

195, 0.8867 0.474 
0.20 120 1901 2.0 0.8862 0.463 
0.25 24 199 act 0.8857 0.443 
0.30 9.2605 50.08 Ip! 0.8856 0.425 
0.35 5.0791 19.98 Zo 0.8859 0.409 
0.40 3.3234 10.44 2.4 0.8865 0.393 
0.45 2.4786 6.46 pie) 0.8873 0.380 
0.50 2.0000 4.47 2.6 0.8882 0.367 
0.55 1.7024 3535 Pag 0.8893 ON355 
0.60 1.5046 2.65 2.8 0.8905 0.344 
0.65 1.3663 2.18 29 0.8917 0.334 
0.70 1.2638 1.85 3 0.8930 0,325 
0.75 1.1906 1.61 Sal 0.8943 0.316 
0.80 1.1330 1.43 Sy) 0.8957 0.307 
0.85 1.0880 129 3.3 0.8970 0.299 
0.90 1.0522 5 i 3.4 0.8984 0.292 
0.95 1.0234 1.08 CYS) 0.8997 0.285 
1.00 1.0000 1.00 3.6 0.9011 0.278 
1.05 0.9803 0.934 Baw 0.9025 0.272 
1.10 0.9649 0.878 Bate 0.9038 0.266 
145 0.9517 0.830 nD 0.9051 0.260 
1.20 0.9407 0.787 “ 0.9064 0.254 
£25 0.9314 0.750 4.10 0.9077 0.249 
1.30 0.9236 0.716 4.2 0.9089 0.244 
135 0.9170 0.687 4.3 0.9102 0.239 
1.40 0.9114 0.660 4.4 0.9144 0.235 
1.45 0.9067 0.635 4.5 0.9126 0.230 
1.50 0.9027 0.613 4.6 0.9137 0.226 
fo) 0.8994 0.593 4.7 0.9149 0.222 
1.60 0.8966 0.574 4.8 0.9160 0.218 
1.65 0.8942 0.556 4.9 0.9171 0.214 
1.70 0.8922 0.540 5 0.9182 0.210 


ie) 0.8906 0.525 aya! 0.9192 0.207 


0.9202 
0.9213 
0.9222 
019232 
0.9241 
0.9251 
0.9260 
0.9269 
0.9277 
0.9286 
0.9294 
0.9302 
0.9310 
0.9316 
0.9325 
0.9333 
0.9340 
0.9347 
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178 Appendices (tables) 
7. LAPLACE TRANSFORMS 


X(t) L(p) 
Dirac impulse (6 function) 1 
1 
Unit step function P 
a 
at ed 
p? 
n! 
t* 
pie 
1 
e7at 
Ise a 
1 
remee —— 
(P +a)? 
‘ a) 
sin wt ————— 
w/p? + w? 
P 
cos wt ne Ste 
Passa 
: 7) 
e-* sin wt eee 
(P + aw)? + 
i Jal (08 
e-* cos wt 
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9. GAMMA LAW 


f(x) = 2° exp (—Ax)X"™ T/n! 
I(n) = (n — 1)! forne N 


foo} 


Po) = | et exp(=tyat 


oc E(@) 

1.0 1.0000 
Lot 0.9514 
12 0.9182 
L3 0.8975 
1.4 0.8873 
iS 0.8862 
1.6 0.893 5 
1.7 0.908 6 
1.8 0.9314 
iF) 0.9618 
2.0 1.0000 
Zl 1.0465 
Dep 1.1018 
2.3 1.1667 
2.4 1.2422 
DS 1.3293 
2.6 1.4296 
Pid 1.5447 
2.8 1.6765 
2.9 1.8274 
3.0 2.0000 
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Index 


ABC (Pareto) analysis 101-4 
Availability definition 10 
determination 38-42 


Bayes theorem 119 
Binomal law 74, 120 


Cause-and-effect analysis 99 ff 
Central limit theorem 127 
Chi-square (x”) test 24, 59, 67, 
108, Table 3 
Confidence interval for mean 
MTBF 130, 131 
Control actions 77 
BASIC program 80 
charts 61-78 
for attributes 73 
individual values 72 
mean 62 
range 69 
standard deviation 67 
levels 90 
Criticality coefficient 47 


Databases 
AVCO 25-32 
CNET 25-32 
NASA 25-32 
RADC 25-32 
US Navy 25-32 
Drifts (in observations) 77 


Experimental design 112-16 
Exponential law 10, 128 


Failure function 
cumulative 10 
instantaneous 10 

Failure Mode Analysis 
(FMA) 45-55 

Fisher test 110, Table 4 

Function 
analysis 47 
tree 52 


Goodness-of-fit tests 24, 59 
chi-square (y’) test 24, Table 3 


Ishikawa diagram 99-101 


Laplace transform 39, Table 7 
log-normal law 12, 127 


Monitoring 61 ff 
Monte Carlo method 42 
Markov chain 36 


. MTBF (Mean Time Between 


Failures) 11 
MTTR (Mean Time To 
Repair) 38 


Ownership costs 9 


Pareto (ABC) analysis 101-4 
Probability laws, distributions 
(definitions, properties) 

binomial 120 
chi-square Table 3 
exponential 128 
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Probability laws (cont. ) 
Fisher Table 4 
hypergeometric 122 
log-normal 127 
normal (Gauss) 124, Table 1 
Poisson 123 
‘Student’ Table 2 
Weibull 129, Table 5 
Probability theory 118 


Quality 
assurance 4 
audit 6 
certification 2, 6 
circles 9 
control 5, 82 ff 
definitions 5 
diagnosis 8 
scales 3 
tools for achieving 7 


Random numbers 43, Table 8 
Rank correlation 104 
BASIC program 105 
Ranks 
method of mean/median 19, 
Table 6 
Redundancy 35 
Regression 
linear 132 
Reliability 
definition, function 10 
laws 10 ff 
models 11 ff 
of systems 32 ff 


Risk analysis 
to customer/supplier 84 


Sampling 


simple, double, multiple 82 ff 


procedures 97-8 
Scatter diagram 116 
Sequential testing 35, 95 
Spearman (rank correlation) 
coefficient 105 
State-change equations, 
graphs 37-42 
‘Student’ test 105, Table 4 


Testing 
sequential 35, 95 
truncated 35 
Transition equations 
graphs 37-42 


Variability 
of products 56 ff 
random 58 
systematic 59 

Variance 
analysis of 106-17 
general linear model 106 
BASIC program 111 
Yates notation 113 


Wald (sequential) test 95 
Weibull law 12-22 
BASIC program 23, 129 


Yates notation 113 
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